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ABSTRACT 

Attention  is  focused  on  the  scalar  reaction-diffusion  equation: 

u  =  A  u  +  f  (u)  ,  where  f  is  cubic-like  and  f(0)  =  f ( 1 )  =  0.  Amongst 
t  x 

spherically  symmetric  solutions  it  is  proved  there  is  a  bounded  unstable 
equilibrium  which  is  decreasing  in  radial  profile.  Under  a  concavity  as. sun;  - 
tion  on  f  this  equilibrium  is  unique.  Moreover  there  is  a  unique  expandin'.: 
spherical  wave  which  is  defined  for  all  time,  positive  and  negative.  As 
t  ->  this  wave  approaches  the  unstable  equilibrium. 

Aronson  and  Weinberger  [2]  have  proved  before  that,  in  all  space  dimen¬ 
sions,  there  are  non-trivial  solutions  that  propagate  (u(x,t)  -*  1  unifornlv 
on  compact  sets  as  t  -*■  +»)  and  ones  that  decay  (u(x,t)  -*■  0  uniformly  as 
t  ->■  +«)  .  This  suggested  the  existence  of  the  unstable  equilibrium. 

There  is  an  interesting  global  description  of  this  propagation/decay  effect. 
The  set  of  initial  data  whose  associated  solutions  approach  the  unstable  equi¬ 
librium  as  t  -*•  +»  splits  a  natural  set  of  functions  into  two  sets.  Data  from 
one  set  yields  a  solution  that  propagates,  and  data  from  the  other  set,  a  solu¬ 
tion  that  decays.  This  fact  is  closely  related  to  the  uniqueness  of  the  expand¬ 
ing  spherical  wave. 
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SIGNIFICANCE  AND  EXPLANATION 


1/ 

Reaction-Diffusion  Equations  have  been  used  to  model  nerve  impulse  pro¬ 
pagation  and  spatially  inhomogeneous  situations  in  chemically  reacting  systems 
and  population  genetics.  The  kind  of  solutions  that  are  often  of  interest  in 
these  areas  are  wave-like  and  do  not  die  out. 

If  the  underlying  spatial  domain  for  the  equation  is  one-dimensional ,  the 
description  of  such  wave  behaviour  amounts  to  finding  a  travelling  wave.  This 
is  a  solution  whose  evolution  in  time  under  the  equation  is  given  by  translating 
along  the  axis.  A  nerve  impulse  is  an  example  of  a  travelling  wave. 

If  a  wave  in  higher  dimensions  is  expected,  for  instance  one  that  is  a 
function  only  of  radius  in  a  spherical  co-ordinate  system,  this  mathematical 
approach  is  not  available.  Such  spherical  waves  look  like  one  dimensional  waves 
a  long  way  out.  In  this  report  we  invert  this  idea  and  deduce  spherical  wave- 
behaviour  in  a  model  example  where  there  is  a  stable  one  dimensional  travelling 
wave.  This  approach  also  determines  the  remaining  work  that  is  necessary  for 
a  complete  picture  of  the  spherical  wave  behaviour  for  this  equation. 
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SPHERICALLY  SYMMETRIC  WAVES  OF  A  REACTION-DIFFUSION  EOLATION 
Christopher  K.  R.  T.  Jones 

Chapter  1 

Introduction 
I .  MOTIVATION 

Reaction  -  Diffusion  equations  often  exhibit  persistent  wave-like 
behaviour.  Examples  of  this  are  the  travelling  pulses  of  the  Fitzhugh- 
Nagumo  and  Hodgkin-Huxley  equations,  travelling  fronts  in  some  scalar 
equations  such  as  the  Fisher  equation  and  target  patterns  and  spiral 
waves  in  some  models  of  the  Belousov-Zhabotinski i  reaction.  Such  phe¬ 
nomena  as  these  usually  depend  on  the  underlying  spatial  domain  being 
unbounded.  This  is  in  sharp  contrast  to  the  case  of  scalar 
reaction  -  diffusion  equations  on  a  bounded  domain  where  energy  argu¬ 
ments  show  that  all  solutions  decay  to  some  equilibrium. 

The  wave  phenomena  in  one  space  dimension,  such  as  the  travelling 
pulses  and  fronts  mentioned  above,  are  often  mathematically  tractable. 
This  approach  is  to  find  an  appropriate  solution  to  the  travelling  wave 
equations,  which  are  ordinary  differential  equations,  and  then  try  to 
prove  that  they  are  stable  as  solutions  to  the  full  partial  differen¬ 
tial  equation.  When  genuinely  higher  dimensional  effects  are  involved, 
as  in  the  case  of  target  patterns  or  spiral  waves,  this  mathematical 
attach  is  not  available.  Results  have  been  obtained  in  some  model 
situations  derived  by  making  simplifying  assumptions,  as  in  the  work  of 
Greenberg,  Hassard  and  Hastings  [1]  and  Kopell  and  Howard  [1,2,31.  The 
main  difficulty  appears  to  be  that,  in  general,  there  is  not  a  specific 
mathematical  object,  like  a  travelling  wave,  whose  discovery  would 
prove  the  existence  and  supply  a  description  of  such  behaviour. 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-75-C-0024 
and  Grant  No.  DAAG29-77-G-0004 . 


In  this  work  we  develop  a  technique  for  handling  spherically  sym¬ 
metric  waves  of  reaction  -  diffusion  equations.  If  a  system  has  s; h  v- 
ically  symmetric  waves  then  it  would  also  have  some  corresponding 
one-dimensional  waves,  loosely  speaking,  these  would  be  the  spherically 
symmetric  waves  "at  infinity".  The  technique  is  to  invert  this  idea 
and  deduce  the  existence  of  spherical  type  wave  behaviour  from  knowing 
one-dimensional  behaviour.  We  hope  that  this  approach  will  give  some 
general  conditions  under  which  such  higher  dimensional  wave  behaviour 
exists  but  it  cannot  give  details  about  the  formation  of  these  waves; 
this  is  an  inherently  higher-dimensional  problem  and  could  not  be 
deduced  from  knowledge  of  the  one-dimensional  mechanisms. 

We  focus  our  attention  on  the  bistable  equation  which  is  described 
in  section  III  below,  since  that  equation  has  the  useful  property  of 
possessing  an  unambigously  stable  travelling  wave.  The  result  of 
applying  the  technique  to  this  equation  is  that  it  guarantees  the  exis¬ 
tence  of  spherical  waves  that  are  asymptotic  to  the  one-dimensional 
travelling  wave  at  infinity  and  also  isolates  the  work  that  needs  to  be 
done  to  get  a  complete  picture  of  the  phenomenon;  we  also  carry  out 
this  work. 

II .  GENERALITIES 

In  the  most  general  case,  we  shall  consider  a  system  of  reaction  - 
diffusion  equations  of  the  form 

ufc  =  Au  +  f (u) 


(1.1) 
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where 


u  e 


x 


n  3^ 

(x  , ...,x  )  e  B  ,  A  =  y  — -  and 

1  n  . u.  „  2 

1=1  3X. 


Au  =  (Au  , ...,Au  ).  We  shall  always  assume  that  f  :  -  lRm  is 

1  m 

smooth. 

The  natural  problem  for  such  a  system  on  all  of  Fn  is  an 
initial-value  problem,  i.e.  a  solution  u(x,t)  should  be  determined  by 
its  initial  condition  u(x) 


(1.2)  u(x,0)  =  u(x) 

We  can  therefore  hope  that  such  a  system  determines  a  semiflow  on 
some  appropriate  function  space.  A  semiflow  on  a  space  Y  is  a  func¬ 
tion  S  :  Y  x  1R+  -*  Y  (whose  domain  may  not  be  all  of  Y  *  ]R+ ,  but 
must  be  an  open  subset)  which  satisfies  (1)  S(S(y,t),s)  =  S(y,t+s) 
and  (2)  S  is  continuous  on  its  domain.  S  is  said  to  be  a  local 
semiflow  if  for  each  y  e  Y,  a  set  of  the  form  (y,[0,s))  for  some 
s  >  0  is  in  the  domain  of  S.  S  is  a  global  semiflow  if  it  is 
defined  on  all  of  Y  x  ]R+. 

For  the  equation  (1.1)  Y  will  be  some  function  space  and 
S(y,t)  =  u(x,t)  where  u(x,t)  is  the  solution  with  initial  condition 
u(x,0)  =  y.  The  space  Y  will  usually  be  the  space  of  bounded  uni¬ 
formly  continuous  functions  u  :  ]Rn  -*•  IR™  which  will  be  denoted  B. 

We  are  interested  in  locating  parts  of  B  that  exhibit  certain  asymp¬ 
totic  behaviour  (in  particular,  wave-like  behaviour)  when  the  semiflow 
is  applied.  The  following  definitions,  due  to  Fife  [1,2],  clarify  this 
quest. 
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The  group  of  transformations  that  leave  the  Laplacian  invariant, 
and  so  also  the  equation,  are  exactly  the  rigid  motions;  call  this 
group  R(n) .  R(n)  acts  in  a  natural  way  on  B;  if  T  c  R(n)  then 
Tu(x)  =  u(Tx) .  Therefore  one  can  see  that  the  space  B/R(n)  inherits 
the  semiflow  which  shall  also  be  denoted  by  S. 

Definition  1.1.  u,v  e  B/R (n)  are  said  to  be  asymptotically  equivalent 
if  there  exists  a  t  e  K  such  that 

lim  ||s(u,t)  -  S(v,t+t)||  =  0 

t-KO 

where  the  norm  used  is  the  sup  norm  over  all  of  Kn  . 

Definition  1.2.  An  asymptotic  state  G  c  B/R(n)  is  an  equivalence 
class  of  initial  data  under  this  relation.  G  is  said  to  be  a  stable 
asymptotic  state  (SAS)  if  it  is  open  in  B/R(n) . 

As  mentioned  earlier  one  attack  on  reaction  -  diffusion  equations 
is  to  look  for  special  solutions  such  as  equilibrium  solutions,  travel¬ 
ling  waves  etc.  In  general  the  least  one  might  expect  of  solutions 
that  warrant  such  attention  is  that  they  be  defined  for  all  t  e  ]R. 
Such  a  solution  would  certainly  be  "special”  as  the  equations  are  para¬ 
bolic  and  so  only  forward  existence  would  be  expected  in  general. 
Definition  1.3.  A  solution  u(x,t)  to  (1.1)  is  said  to  be  a  permanent 
solution  if  it  is  defined  for  all  t  e  K. 

Satisfaction  of  the  ultimate  dream  for  a  given  reaction  - 
diffusion  equation  would  be  a  description  of  all  its  SAS's  and  the 
location  of  all  their  bounded  permanent  solutions. 

We  cannot  usually  expect  to  be  able  to  prove  that  a  certain  set  is 
a  stable  asymptotic  state  with  respect  to  all  of  B,  but  in  a  given 


problem  there  may  be  a  natural  subspace  of  B  in  which  the  asymptotic 
state  can  easily  be  proved  to  be  open.  This  will,  only  be  meaningful  if 
that  subspace  is  invariant,  or  at  least  positively  invariant. 

Definition  1.4.  A  c  y  is  invariant  under  S,  if  S(t)A  =  A  for  all 
t  >_  0.  A  is  positively  invariant  if  S(t)A  c  a  for  all  t  21  0  ■ 

Suppose  A  c  b  is  positively  invariant. 

Definition  1.5.  A  set  G  is  an  SAS  with  respect  to  A  if  it  is  an 
asymptotic  state  and  it  is  open  in  A. 

In  regard  to  looking  for  permanent  solutions  one  of  the  useful 
ideas  of  dynamical  systems  is  that  of  an  co-limit  set.  If  A  c  y  then 

W  (A)  =  H  cl  (A  •  [t,°°)) 
t>0 

The  crucial  property  of  w-limit  sets  is  that  they  are  closed  and 
invariant. 

A  consequence  of  this  is  that  points  in  io(A)  have  local  backward 
existence  but  the  solutions  may  blow  up  in  finite  backward  time,  so 
to { A)  need  not  consist  entirely  of  permanent  solutions.  If  however  we 
knew  that  S  was  a  compact  semiflow  on  Y  (see  chapter  4,  section  I) 
then  to ( A)  is  a  compact,  non-empty  set  that  consists  entirely  of  per¬ 

manent  solutions. 

This  compactness  property  is  not  satisfied  on  B  if  it  carries 
the  sup-norm  topology,  however  we  can  replace  this  with  the  compact- 
open  topology,  without  destroying  the  semiflow  property,  and  the 

* 

for  notation  see  section  IV 


compactness  condition  is  satisfied.  Permanent  solutions  s', 
found  by  calculating  the  u-liir.it  set  of  a  bounded  connect- - 
the  SAS  in  question. 


III.  THE  BISTABLE  EQUATION 

Although  the  main  abstract  construction  of  this  work  a 
any  system  of  reaction  -  diffusion  equations  the  only  appli 
worked  out  is  tc  the  scalar  equation/  called  the  bistable  e 

(1.3)  u  =  Au  +  f(u)  / 


where  u  e  ]R  .  We  split  the  assumptions  on  the  nonlinearit 


two  parts: 


(HI) 


f  :  R  -*•  IR  is  smooth, 
unique  a,  between  0 
Furthermore  f ' (0)  <  0 


f(0)  =  f (1)  =  0  and  the 
and  1  such  that  f(a) 

f 

f ' (1)  <  0  and  /  f (u)du 


0 

(H2)  For  all  a  <  B  <  1  f"(S)  <_  0. 

Roughly  speaking,  if  f  satisfies  (HI)  then  it  is  of 


Figure  1.1 


where  the  area  of  the  hump  above  the  x-axis  is  biager  than  that  of 
hump  below.  If  f  satisfies  (H2)  then  the  hump  above  is  concave  down. 

It  will  always  be  assumed  that  f  satisfies  (HI),  (H2)  will  be  use':  m 
a  crucial  uniqueness  result  in  chapter  2  and  it  will  be  clear  w.at 
elusions  can  be  made  when  (H2)  is  not  satisfied.  We  shall  assume 
nothing  about  f  outside  the  interval  [0,1]  as  all  our  attention 
will  be  restricted  to  data  in  this  set;  this  is  possible  as,  by  the 
maximum  principle  (see  Chapter  3,  section  II):  if  0  £  u(x,0)  <_  1 
then  0  <_  u(x,t)  £  1  for  all  t  £  0. 

The  crucial  property  of  the  scalar  equation  with  this  f(u)  is 
that  it  has  a  unique  (up  to  translation)  travelling  wave  in  one  space 
dimension  which  is  stable  in  a  very  strong  sense.  In  other  words,  it 
has  a  unique  solution  of  the  form  u(x-ct)  which  satisfies  the  bound¬ 
ary  conditions  u(-“)  =  1  and  u(+»)  =  0  (see  Aronson  and  Weinberger 
[1]  and  Fife  and  McLeod  [1]).  Fife  and  McLeod’s  theorem  is  the  stron¬ 
gest  stability  statement,  it  says  that  the  travelling  wave  is  exponen- 
tailly  stable  to  all  initial  data  which  lie  above  the  middle  zero  at 
“°°  and  below  it  at  +°°  (it  is  reproduced  here  as  our  theorem  5.1)  . 

A  standard  maximum  principle  argument  implies  that  u  =  0  and 
u  i  1  are  stable  (this  is  the  source  of  its  name)  to  perturbations  in 
the  sup-norm,  u  =  a  is  unstable  by  the  same  reasoning.  The  travel¬ 
ling  wave  determines  the  asymptotic  behaviour  of  data  that  are  transi¬ 
tions  between  the  two  stable  states.  Heuristically  we  could  take  for 
initial  data  characteristic  functions  of  a  set;  then  if  u(x,0)  is  the 
characteristic  function  of  a  half-line  (-°°,a] ,  u(x,t)  tends  to  a 
translate  of  the  travelling  wave  as  t  -*•  +°°. 

ja 


In  any  space  dimension  we  could  take  u(x,0)  to  be  the  character¬ 
istic  function  of  a  compact  set  D.  If  D  is  large  enough  and  more  or 
less  spherical  then  at  the  edge  of  D  u(x,0)  will  not  be  too  unlike  a 
plane  front  and  so  one  might  expect  it  to  propagate  in  the  same  fashion 
as  a  one-dimensional  wave.  On  the  other  hand  if  D  were  too  small  we 
might  expect  the  solution  to  tend  to  zero  everywhere.  Aronson  and 
Weinberger  (2]  have  shown  that  both  of  these  phenomena  occur.  They 
give  an  integral  condition/  that  will  cover  some  data  of  the  above 
form  if  D  is  small  enough,  which  ensures  decay  of  the  corresponding 
solution.  They  also  give  a  condition  for  propagation  (u(x,t)  -*■  1 
uniformly  on  compact  sets) .  The  latter  theorem  specialized  to  the  case 
under  consideration  is: 

Theorem  1 . 1  (Aronson  and  Weinberger  [2]).  With  f  satisfying  (Hi) 
there  is  an  indexed  family  of  functions  on  ]Rn  say  (u^|i  c  A},  each 
with  compact  support,  so  that  if,  for  some  \ 

u(x,0)  >  u.  (x)  for  all  x  e  lRn 
A 

and  u(x,t)  is  a  solution  of  (1.3)  then  u(x,t)  ->  1  uniformly  on  com¬ 
pact  sets. 

They  also  prove  that  the  propagation  has  an  asymptotic  speed  which 
is  that  of  the  one-dimensional  travelling  wave,  but  we  shall  not  go 
into  this  here. 

In  particular,  from  the  above  theorem  we  know  that  there  is  spher¬ 
ical  wave  propagation  for  (1.3),  it  is  this  pehnomenon  that  we  shall 


analyze . 


XV.  SYNOPSIS 
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From  Aronson  and  Weinberger's  results,  there  is  a  'threshold 
effect',  i.e.  some  solutions  decay  and  some  propagate;  this  suggests 
that  there  is  an  unstable  equilibrium  solution  which  demarcates  the 
boundary  between  these  two  regimes  of  behaviour.  In  chapter  2  we  prove 
the  existence  of  such  solutions  and  under  the  assumption  that  f 
satisfies  (H2)  we  prove  it  is  unique.  This  solution  is  obviously 
unstable  and  under  (H2)  we  prove  that  the  linearized  operator  at  this 
equilibrium  solution  has  one  positive  eigenvalue  which,  loosely  speak¬ 
ing,  represents  one  dimension  of  decay  and  propagation. 

The  remainder  of  the  work  is  devoted  to  giving  a  global  picture  of 
this  propagation/decay  behaviour,  using  the  solution  found  in  chapter  2 
as  a  pivot. 

In  chapter  3  we  compile  the  necessary  information  for  applying  the 
concepts  and  methods  of  dynamical  systems.  Except  for  section  IV  this 
chapter  is  all  standard  material. 

In  chapter  4  we  construct  the  basic  machine  that  gives  out  spheri¬ 
cal  information  when  we  plug  in  facts  about  the  one-dimensional  behav¬ 
iour.  In  this  chapter  we  describe  the  "spherical  attractor"  when  f 
satisfies  (H2) ,  this  is  the  set  of  permanent  solutions  that  are  rele¬ 
vant  to  this  propagation/decay  effect. 

In  chapter  5  we  refine  the  machine  of  chapter  4  to  include  a  mov¬ 
ing  co-ordinate  frame;  this  gives  finer  information  about  the  spherical 
waves.  We  also  prove  the  necessary  one-dimensional  facts  in  this 
chapter . 


In  chapter  6  we  sum  the  work  up  and  translate  the  r<  .  si*:  r :  - 
language  of  dynamical  systems  into  terms  that  are  cl  >«vi  •  • 

ture  on  reaction  -  diffusion  equations. 

Notation 

(1)  For  a  set  A:  Int(A)  =  interior  of  A; 

cl (A)  =  closure  of  A;  A®  =  complement  of  A. 

(2)  If  S  :  Y  x  ]R+  -*■  Y  is  a  semiflow  and  A  c  Y,  I  r  1R+ 

S (t) A  =  { S ( t ) y | y  €  A};  A  •  t  =  S (t) A; 

A  *  I  =  {y  •  t | y  e  A,  t  c  I}. 

(3)  B  =  {u  :  IRn  -*■  IRn  |u  is  bounded  and  uniformly  continuous- 

HUIIB  =  HUIL  =  Sup  lu(x)  I 

_n  *  •••*•* 

►  •  •  *  *  xtlR 

n 

(4)  If  P  =  TT  [a.,b.]  is  an  invariant  rectangle  (see  chapter  3, 

i=l  1  X 

section  2) ,  then  M(P)  c  b  is  given  by: 

M(P)  =  (u(x)  =  (u(x)  =  (u  (x) , . . . , u  (x))|a.  <  u . (x)  <  b.  : 

1  m  i—i  —  i 

every  x  e  lRn}  and  supplied  with  the  compact-open  topoloq 
The  dependence  on  P  is  usually  omitted. 

(5)  If  u(x,t)  and  u(x)  are  used  together,  then  u(x,0)  =  u(x). 


I.  ONE  DIMENSIONAL  CASE 

In  this  chapter  we  shall  prove  theorems  on  the  existence,  unique 
ness  and  spectral  properties  'of  nonconstant  equilibrium  solutions  of 
(1.3).  We  shall  firstly  review  the  well  known  case  of  one  s:  ace  dime 
sion. 

An  equilibrium  solution  in  one  space  dimension  satisfies  the 
equation 

(2.1)  u  +  f  ( u)  =  0  . 

XX 

We  are  interested  in  solutions  that  decay  to  0  '  at’  x  =  r ■ .  Trans- 
forminq  (2.1)  into  a  system  of  ODE's  gives 


(2.2) 


u'  =  v 

•  =  A. 

dx 

v'  =  -f(u) 


which  is  Hamiltonian  and  so  the  phase  portrait  is  easily  drawn 


11 


*  *  •  Figure  2.1 


From  this  we  can  see  that  there  is  a  unique  non-trivial  solution  (up  to 
translation)  (u(x),v(x))  which  satisfies  the  boundary  conditions, 
namely  the  "body  of  the  fish". 

Looking  for  the  eigenvalues  of  the  linearization  around  this  solu¬ 
tion  leads  to  the  equation 

u '  =  v 

<2.3>  1  -  S  • 

v'  =  (X  -  f'(u))u 

Using  the  fact  that  (v(x),v'(x))  is  a  solution  for  1=0,  by  stan¬ 
dard  comparison  arguments  it  can  be  shown  that  there  is  a  unique 
A  >  0  for  which  (2.3)  admits  a  bounded  solution. 


In  the  following  sections  we  shall  show  that  the  same  structure 
shows  up  in  the  higher  dimensional  case. 


II.  EXISTENCE 


We  seek  spherically  symmetric  equilibrium  solutions  of  (1.3); 
these  are  solutions  of 


u  +  -  u  +  f (u)  =  0 

rr  r  r 


Theorem  2.1.  For  fixed  n  and  f  satisfying  (HI)  there  is  at  least 
one  bounded  non-constant  solution  of  (2.4)  that  satisfies 


u(r  =  +“)  =0 


Proof ;  Transforming  (2.4)  into  a  system,  we  get 


v'  =  -—-v  -  f(u)  . 


Roughly  speaking,  when  r  =  +°°  we  have  the  equilibrium  equation 

for  the  one-dimensional  problem.  In  order  to  make  this  precise  we  per- 

r 

form  the  transformation  p  =  — — 


(2.6)  V  =  -  _  f(u)  •  =  A. 

p  -  dr 

p'  =  (1-P)2  . 


There  is  a  singular  surface  at  p  =  0,  but  otherwise  this  equa¬ 


tion  is  well  defined  on  R  *  [0,11  .  The  phase  portrait  in  the 


section  p  =  1  is  easily  seen  to  be  that  of  Figure  2.1. 


Since  the  solution  must  be  regular  at  the  origin,  w< 


tion  of  (2.6)  that  satisfies  the  boundary  conditions 


(2.7) 

(2.8) 


v(0)  =  0 
u(°°)  =  0 


(2.8)  means  that  the  purported  solution  must  lie  in  the  set  W  civ-:, 
by: 

W  =  {(u,v,p)  |  (u(r) ,v(r) ,p (r) )  -*  (0,0,1)  as  r  -*•  +* 

and  (u(0) ,v(0) ,p (0) )  =  (u,v,p) } 

By  standard  results,  see  for  instance  Kelley  [1],  there  is  a  1 

c  s 

local  center-stable  manifold,  call  it  W,  ,  at  (0,0,1)  whose  tan- 

loc 

gent  space  is  generated  by  £  and  (0,0,-l)  where  £  is  anv  vector 

s  s 

.  ..  .  *  ....  •  •  •  *  ‘ 
in  u , v  space  tangent  to  the  stable  manifold  of  (0,0)  in  Figure  2.1. 

Now  let  T(r)  be  the  solution  operator  of  (2.6),  for  all  r 

T(r)W<Ts  c  w.  Define  WCS  =  U  T(r)W^s  ,  then  WCS  c  w.  Since 
loc  _  loc 

r<0 

cs  _ 

W  n  {p=l}  contains  points  in  the  positive  quadrant  of  the  u,v 


plane,  see  Figure  2.1,  if  r  is  sufficiently  large  negative  T(r)W 


loc 


,cs 


contains  these  points.  Therefore,  since  W  "  is  a  manifold  and  T(r) 

loc 

CS  - 

is  a  dif feomorphism,  T(r)W,  n  {p  =  p  <  1}  intersects  the  positive 

loc 

quadrant  of  p  =  p,  for  p  close  to  1. 

The  points  W  n  {u,  v  >  0}  stay  in  {u,  v  >  0}  under  application 
of  T(r)  for  r  <  0.  To  see  this  we  must  consider  the  Hamiltonian  of 
the  one-dimensional  equation 

v2  .u 

H (u, v)  =  —  +  /  f (s) ds  . 

0 


i 


We  compute  H  along  orbits  of  (2.6) 

A  =  -lnzlMlzoIv2 
P 

so  if  p  ?  1  and  v  ^  0  then  A  <  0.  For  the  one-dimensional  equa¬ 
tion  the  energy  curves  are  the  solution  curves.  If  p  /  1,  the  solu¬ 
tions  of  (2.6)  cross  these  curves  with  decreasing  energy.  We  borrow 
the  two  curves  and  in  Figure  2.2  belcw  from  Figure  2.1  (they 

are  part  of  the  energy  surface  H  =  0)  and  setting 

=  {v=0,  y  <_  u  <_  1}  and  =  {u=l,  v  >  0},  it  is  easy  to  see  that 
the  vector  field  in  the  p  =  p  slice  is  as  depicted  on 

C  =  C  t  C_  U  C  U  C .  . 

12  3  4 


Figure  2 . 2 


so  the  region  bounded  by  C  in  the  upper  half-plane  is  negatively 
•  cs 

invariant.  Since  H  <  0,  W  must  lie  outside  H  =  0  and  so  it 
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crosses  C  on  for  p  close  to  1.  Therefore  WCS  intersects 

the  positive  quadrant  for  every  p  >  0. 

Let  e  >  0  be  a  fixed  snail  number  and  define  a  set  D  bounded 
by  the  four  lines 

u  =  a,  -z  <_  v  <_  e;  v  =  a  ’ <_  u  £  1 ;  u  =  1 ,  -e  <_  v  <_  e 

From  the  above  W”J  o  D  r.  {p=p}  ^  $  for  every  p  >  0.  We  want  to  show 

cs  — 

there  is  a  solution  in  W  which  is  in  D  at  some  p  and  stays  in 
for  all  0  <_  p  <_  p ,  renormalizing  the  independent  variable  then  gives 
one  that  satisfies  (2.7). 

Let  D  =  3D  n  { v>0}  and  D  =  3D  n  {v<0},  then  for  every  o, 

wc~  n  {p=p}  n  D  *  $  and  wcs  n  {p=p}  n  D  /  tp. 

U  L 

If  p  is  small  enough  the  vector  field  on  3D  is  as  shown  below 


Figure  2.3 


For  this  p,  let  F  be  a  compact  curve  in  WCS  n  {p=p}  that  inter¬ 
sects  Dy  and  Dy,  such  a  curve  clearly  exists.  Let  Fy  =  {y  e  F  j 

for  some  r  <  0  T(r)y  4  D  but  T(r)y  e  {v>0}}  and  F  be  the 

L 
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corresponding  set  with  v  <  0.  Fy  and  F^  are  obviously  both  open 
and  from  the  picture  in  Figure  2.3  are  both  non-empty  but  then 

F  *  Fu  11  fl 

so  there  exists  y  t  F  that  stays  in  B  for  all  0  such  that 
0  <_  p  d  ■  This  completes  the  proof  of  Theorem  2.1. 


This  solution,  u(r),  found  in  theorem  1  satisfies  the  following 
two  properties 

(A)  u'(r)  <  0  for  all  r  e  (0,+°°). 

(B)  y  <  u(0)  <1  where  Y  :  a  <  Y  <  1  is  determined  by  the 

Y 

condition  /  f(s)ds  =  0. 

0 

Property  (A)  is  clear  from  the  fact  that  the  region  described  in 

cs 

Figure  2.2  is  negatively  invariant.  Since  W  lies  outside  H  =  0 
it  is  clear  that  u(0)  >  y.  To  see  that  u(0)  <  1  we  apply  the  maxi¬ 
mum  principle.  Since  u  =  1  satisfies  (2.4),  by  an  application  of  the 
mean  value  theorem  w  =  1-u  satisfies 

(2.9)  Aw  +  c(x)w  =  0 

for  some  bounded  function  c(x).  But  by  (A)  and  the  obvious  fact  that 
u(0)  _<  1  we  see  that  1-u  _>  0;  but  then,  by  the  maximum  principle, 
see  Protter  and  Weinberger  [1] ,  1-u  cannot  have  an  interior  minimum 

of  zero  which  it  would  if  u(0)  =  1, 


III.  SPECTRAL  ANALYSIS 


In  this  section  we  consider  the  eigenvalue  p  rob  lor.  for  th-  ;.u  . 

rized  operator  around  a  solution  u  of  the  form  found  ir.  :••••• 
Specifically  we  want  to  know  for  which  A's  the  equation 

(2.10)  Au  +  (f 1 (u)-\)u  =  0 

admits  a  bounded  solution.  In  particular  we  want  to  know  for  which 
A  >_  0  (2.10)  admits  a  bounded  solution. 

It  follows  from  standard  results,  see  e.g.  Reed  and  Simon  [1]  v. 
IV,  that  any  positive  eigenvalue  must  correspond  to  a  spherically  sy 
metric  eigenfunction.  This  comes  from  the  fact  that  o  is  an  eigen¬ 
value  with  an  associated  n-dimensional  eigenspace  generated  by 
{u  ,...,u  }  (this  is  the  translation  eigenvalue)  and  each  of  these 

X1  Xn  ... 

..  .  H  .  .  ••  ... 

has  no  zero  except  at  zero.  The  rest  of  the  spectrum  is  bounded  away 
from  zero  to  its  left. 

So  we  look,  for  solutions  of  the  equation,  with  A  >_  0, 

(2.11)  u  +  —  u  +  (f '(u)  -  \)u  =  0  . 

rr  r  r 


Theorem  2.2.  If  f  satisfies  (HI)  and  (H2) ,  there  is  only  one  >  ^ 
for  which  (2.11)  admits  a  bounded  solution  and  that  A  is  positive. 
Proof:  Converting  (2.11)  to  a  system 


(2.12) 


u '  =  v 

v'  =  - ~~ v  +  (A  -  f'(u))u  . 


The  associated  angular  equation  is 


0 
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where 


0  =  arctan{— ) 
u 


We  only  have  to  worry  about  solutions  of  (2.13)  . 


Roughly  speaking,  at  r  =  +°°  (2.13)  is  the  angular  equation  for 


the  system 


(2.14) 


v'  =  (A-f ' (u) ) u 


This  equation  has  stable  and  unstable  subspaces  which  are  lines  of 


slope  --'A-  f 1  ( 0 )  and  +v'A  -  f '  ( 0)  respectively.  Let  6  be  the 
angle  (in  the  fourth  quadrant)  that  this  stable  subspace  makes  with  the 

u-axis.  It  is  not  hard  to  see  that  if  (2.12)  has  a  bounded  solution, 

*  *  *  *  *  *  *  •  •  ••  • 

the  corresponding  solution  9(r)  to  (2.13)  must  satisfy 


(2.15) 


V)  =  9  A 


(mod  it) 


It  is  also  standard  that  there  is  a  unique  solution,  0^  (r)  (mod  tt)  , 
of  (2.13)  that  satisfies  (2.15). 

The  proof  requires  a  study  of  0.  (r)  and  its  limit  as  r  tends 
to  0.  The  solution  of  the  full  equation  remains  bounded  only  if 
0.(0)  =  r'  (mod  tt )  ,  so  the  question  is  for  what  A's  can  we  have 


(2.16) 


9  (0)  =  itur  for  some  m  ? 


We  will  split  the  proof  up  into  three  parts  but  firstly  we  make 


some  preliminary  observations.  We  must  have  m  ^  0  as  when  0  =  -tt/2, 


0'  =  -1.  The  A's  which  work  for  m  are  decreasing  in  m,  that  is. 


if  A,,  works  for  nu  i  *  1,2  and  m^  <  m2  then  A 2  <  A  .  This 


follows  from  the  fact  that  if  X  >  A,  and  0^(r),  6.-(r)  are  solutions 
of  (2.13)  with  X  and  A  respectively  then  if  0^(R)  l_  5{(R)  (we  may 

even  have  R  =  +°°  here)  then  9  (r)  <  6-(r)  for  all  r  £  R.  This 
latter  fact  is  proved  by  a  comparison  between  the  X  and  A  eauations. 

From  this  we  know  that  the  set  of  A's  qivinq  an  affirmative 
answer  to  (2.16)  are  bounded  above  and  the  associated  m's  are  increas¬ 
ing  with  decreasing  A.  The  rest  of  the  proof  is  split  up  into  proving 
the  following  three  parts. 

(1)  For  a  given  m  there  is  at  most  one  A  which  satisfies 
(2.16)  . 

(2)  There  is  a  A  >  0  for  which  6.(0)  =  0. 

A 

(3)  0  <  0  (0)  <  tt. 

Once  these  are  proved  the  proof  is  complete  since,  by  (3)  and  the 
decreasing  property,  if  (2.16)  can  be  solved  for  m  =  1,  the  associated 
A  must  be  negative  and  therefore  so  must  all  other  A's  be  negative 
associated  to  any  m  ^  1.  The  only  possible  A  >_  0  must  therefore 
correspond  to  m  =  0,  from  (1)  there  is  a  unique  such  A  and  (2)  shows 
that  it  exists.  The  hypothesis  (H2)  is  used  in  the  proof  of  (3)  and 
this  is  perhaps  the  heart  of  the  theorem. 

Proof  of  (1) :  Suppose  A  satisfies  (2.16)  for  some  m,  we  firstly 
show  that  0^(r)  is  the  only  solution  (mod  7i)  of  (2.13)  (with  A 
fixed  at  this  value)  which  satisfies  (2.16). 

Suppose  it  were  not,  then  we  would  have  two  solutions  (u^,v^) 
and  (u2,v2)  of  (2.12)  with  the  property  that  v^  -  v2  -*■  0  as 
r  -+  0,  and  also  that  u  =  u  -  u2,  v  =  v  -  v2  is  a  solution  of 


21 


(2.12).  The  angular  and  radial  equations  of  (2.12)  are  given  by  (2.13) 

and 

(2.17)  p<  =  -p (sin  9  cos  9(1  +  X  -  f'(u))  -  — sin'  9) 


where 


Figure  2.4 


If  a  solution  remained  in  this  cone  as  r  -*■  0,  from  (2.17)  it  would 
have  either  pHO  or  p-»+«  as  r  -»■  0.  But  by  assumption  we  could 
choose  (u^v^  (u2'v2}  50  that  tu,v>  iies  in  this  cone  for 

some  small  r.  So  if  *  C  we  would  have  to  have  p  =  0  and 

SO  (Uj.»V  )  *  (u2'v2)* 


Assume  that  X  <  •  then  there  must  be  an  k 

6  (R)  >  0  (R)  .  3ut  then  there  exists  a  solution  t  .. 
A1  '2 

>.  =  X  ,  call  it  6.  (r)  ,  sc  that 

2 


9,  (R)  "  v,  (R)  •  •  (R) 

1  ‘2  '2 


but  then  for  all  r  <_  R  e  (r)  >  9  (r)  and  since  1  : 

1  *2 

unique  solution  of  (2.13)  with  X  =  >  that  satisfies 

0.  (0)  >  mit  and  so  6.  (0)  >  mt. 

A2  1 

Proof  of  (2)  :  That  such  a  >,  exists  for  m  =  0  follows  : 
straightforward  shooting  argument  using  the  two  facts:  (i 
then  6A(0)  £  0X  and  (ii)  if  X  =  0,  8^(0)  >  "/2. 

(i)  follows  from  the  fact  that  if  '>>0,  8  =  ■?,  im;  i 
(from  (2.13))  and  so  @A(r)  £  6  for  all  r  1  0.  For  (ii) 
be  compared,  when  X  =  0,  to  the  angular  equation  of 

u'  =  v 


(2.18) 


.  n-1 

v'  —  v 

r 


f'  (u)u  + 


which  has  a  known  solution,  namely  u'(r).  It  is  easy  to  see  that  ts 
forces  0  (0)  ir/2. 

Proof  of  (3) :  For  the  remainder  of  this  section  we  shall  be  discuss i> 
(2.13)  with  X  =  0  and  so,  for  the  sake  of  notation,  will  drop  the 
subscript  0  on  9^  and  0^(r).  With  X  =  0,  (2.13)  is 

(2.19)  0'  =  -  —  —  sin  0  cos  9  -  f‘  (u)cos^  0  -  sin^  9 

r 


s 


-  ...  ■  ■  ■  _ _ 


This  is  also  the  angular  equation  of  the  tangent  (or  first  variar icr.i’ 
equation  of  (2.12)  about  the  solution  (u(r),v(r?)  and  0  is  the 
angle  of  the  stable  nanifold  to  (0,0)  in  Figure  2.1.  Consequently 


6(r)  is  the  angle  cf  the  tangent  line  to  the  manifold 
R  , 


VJ  =  W  n  (p  =  zrr)  at  (u  ( K)  ,  v  ( R)  )  .  v.'e  actually  have  r.nt  shown  that 
R  Ml 


this  is  a  manifold  yet,  a  proof  cf  this  is  in  section  III,  hut  wo  a; 


know  it  contains  a  .manifold ,  namely  W  n  t?  =  — 


K+l  * 


»nr.c*  fcr  tr*c 


following  arguments  this  suffices. 


We  must  prove  that  0(0}  <  n,  we  do  this  in  two  steps,  it  is  in 
step  2  that  the  hypothesis  (H2)  is  used. 

Remark:  (3)  is  about  8(r)  which  is  a  solution  of  (2.13)  asso¬ 


ciated  to  a  u(r)  which  satisfies  (2.4),  (2.8)  and  its  derivative 
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So  then  there  are  points  (u  ,v  )  and  (u  ,v  )  on  W  so  that  if 

XX  2  2  R 

(ik  (r)  ,v^  (r) )  ,  i  =  1,2  satisfy  (ir(R),v  (R))  =  (ir.vj,  we  have 


H(U1'V1>  ?  h^U2'V2) 

and  yet  (u^  (r)  ,v^(r)  )  -*  0  as  r  -»  +°°.  Recall  that 

(n-1)  (v.  (r))2 
ft(u  (r)  ,vi(r)  )  - - — - 


fron  the  definition  of  R  and  the  fact  that  u^  <  we  have 
v^fr)  <  v2<r)  and  so  H(u^(r) ,v^(r) )  <  H (u2 (r) ,v2 (r) )  but  then 
H  (u^  (°°)  ,  v^  (°°)  )  >  H  (u2  (»)  ,  v2  (°°) )  which  is  a  contradiction  and  so 
u(R)  >_  a. 

Step  2:  If  9  (R)  =  7T  then  v(R)  >0  (R  >_  0)  . 

Consider  the  angular  variation  of  the  full  nonlinear  equation 
(2,5)  about  the  point  (a,0)  in  the  u,v-plane, 


ijj  =  arctan(- — ) 
u-a 


This  quantity  satisfies 

(2.20)  <|i'  =  -^y^sin  cos  -  f '  (u)cos2  -  sin2  <J>  +  g(r) 


where  g(r)  is  given  by  the  expression 

f (u) (u-a) 


g(r)  =  -  ”2  +  f'  (u)cos‘  1 |< 

v  +(u-a) 


f (u) (u-a)  +  f 1  (u) (u-a) ' 

-2-2 
v  +(u-a) 


u-a 


-2  -  2 
v  +(u-a) 


(f '  (u)  (u-a)  -  f  (u) ) 


and  so  1  f 


Now  if  f "  (u)  <_  0  then  f'  (u)  (u-a)  -  f(u)  C 
have 

(2.21)  g(r)  <_  0  . 

Suppose  0  ( Rg )  =0,  then  u(R^)  L  ^  and  so  !  •'  '■  : 

(unless  v(r)  becomes  fairly  positive,  which  case  does  r.ot  in*-  :  • 

us).  By  (H2)  if  u  a,  f"(u)  <_  0  and  so  if  r  <_  R  (2.21''  1:  . 

Since  u(RQ)  2.  a»  *MRg)  —  -7r/2- 

0(r)  -  7T  satisfies  (2.19)  and  'ji(r)  satisfies  (2.20)  and  sir..-. 

0<R  )  -  tt  <_  'J'(Rg)  we  must  have  6  (r)  -  n  <_  \p(r)  for  all  r  sue';. 

0  i  r  £  Rq.  So  if  0  (R)  =  tt,  ip  (R)  _>  0  and  so  v(R^)  _>  0 .  If  f’^  * 

we  have  strict  inequality  and  since  0(RQ)  -  n  <  f (R  ) ,  as  bvfot" 

(in  the  proof  of  (1))  i|i(r)  can  be  compared  to  another  solution  of 
(2.19)  and  so  vp ( 0)  >0  in  this  case  also.  In  either  case  v(R)  ■ 
and  the  proof  of  step  2  is  complete. 

To  complete  the  proof  of  (3)  notice  that  we  have  shown  that  if 

0  (R)  =  ir  then  v(R)  >  0,  but  this  is  impossible  for  any  solution  of 

theorem  2.1  as  according  to  property  (A)  v(r)  <_  0.  So  clearly  we  must 
always  have  0(0)  <  tt. 

IV.  UNIQUENESS 

From  the  fact  that  the  region  bounded  by  C  in  Figure  2.2  is 
negatively  invariant  the  only  solutions  to  (2.5)  that  satisfy  the 
boundary  condition  u(+°°)  =  0  and  stay  between  0  and  1  must 
satisfy  properties  (A)  and  (B)  at  the  end  of  section  I. 


Mlt'.,  1JP 


then  let  C  =  y  +  C_.  Let  <f(r)  be  the  solution  of  (2.22)  such  that 
y  0 

y  =  <fi(0)  ,  i|)(r)  be  that  with  z  =  tp (0)  . 

Lemma  2.1.  There  exists  a  neighbourhood  U  of  (0,0,0)  so  that 

if  z,y  £  U  and  z  £  then  |ir^(^(r)  -  )  |  is  increasing  so 

long  as  ip  (r)  and  v’(r)  are  both  in  U. 

Notice  that  (0,0,1)  of  (2.5)  is  transformed  into  (0,0,0)  of 

cs  cs 

(2.22).  For  the  sake  of  notation  call  the  image  of  W,  also  W, 

loc  loc 

I  CS  cs 

Take  z  £  U  but  z  j  W,  ,  clearly  there  is  a  y  £  U  n  W,  so  that 

T  loc  loc 

z  £  Cy,  but  since  <p(r)  -*■  0  as  r  ■*  +°°,  we  cannot  have  i|i(r)  -*-0 


.13 


without  it  leaving  U  first,  by  the  lemma.  Transforming  this  state- 

cs  ,  .  . 

ment  back  to  the  original  variables  gives  us  that  W  =  W  as  desired. 
Proof  of  Lemma  2.1. 

Take  U  to  be  a  ball  of  radius  e  >  0,  so  small  that 
|g(w)  -  g  (w)  |  <_  fi|w-w| 

for  w,w  e  U.  If  P(r)  is  the  evolution  operator  of  (2.22)  then  5  can 
be  chosen  so  that 


P (r)  (C  n  U)  c  C  .  , 
y  'r(r) 

for  y  c  U  so  long  as  V>(r)  e  U;  i.e.  these  cones  are  invariant  under 
the  flow.  Suppose  z  e  then  z-y  e  CQ  and  the  associated  solution 

ip(r',  -  <p(r)  satisfies  the  equation 


(2.23)  =  AOl-t)  +  g W  -  g(<rt 

and  since  A^  >  0,  it  is  clear  that  if  6  is  small  the  function  on 
the  right  points  inward  on  3CQ  n  U. 

The  solution  ^(r)  is  given  by  the  variation  of  constants  formula 

r 

.  ,  Ar  f  A (r-s)  ,  ,  .  ,  . 

Y?(r)  =  e  y  +  J  e  g(^(s))ds 

0 


and  so, 


A  r  r  A  (r-s) 

it  (tp(r)  -  ss(r))  =  e  tt  (z-y)  +  |  e  n  (g(ij)(s))  -  g(v>(s))) 

0 


ds 


Because  tp  ( s )  -  y?( s)  e  CQ  there  must  exist  a  k^  >  0  so  that 


llr1(g(<l'(s))  -g(V>(s)))|  _<  k1 6|m1(4i{s)  -  ^(s) ) 


9 


but  then 


X  r  r  X  (r-s) 

|  tt  (iMr)  -  s5(r))|  >•  e  [  rr  (z-y)  |  -  /  e  k  6 1  t  (s)  -  y'(s)  )  !  ds 

1  0  1  1 

Now  if  we  assume  that  |t1(^(s)  -  s?(s)  )  |  <_  k2e>,S  j  it  (z-y)  |  where 
X  <  (otherwise  we  would  be  done)  we  get 

X.r 

|  tt^  ( ^  ( r )  -  s^(  r ) )  |  _>  e  (|n  (z-y)|  -  k6|ir  (z-y)|) 

V  .  , 

=  e  ( l-k6) | n  (z-y) | 


If  6  is  chosen  small  enough  the  lemma  follows. 


We  can  now  state  and  prove  the  uniqueness  theorem. 

Theorem  3.3.  If  f  satisfies  (HI)  and  (H2)  then  the  nonconstant  solu¬ 
tion  of  Theorem  3.1  is  the  unique  one  between  0  and  1. 

Proof.  We  firstly  fix  some  notation.  From  lemma  3.1  we  know  that  for 


each  R,  W  is  a  one-dimensional  curve,  and  so  W  \(u,v)  has  two 
R  R 


connected  components  call  these  W^(u,v)  and  W^(u,v),  where  they  are 

R  R 


determined  by  the  requirement  0  €  W  (u,v) .  Let  T  (u,v)  be  the  tan- 

R  R 


gent  line  to  W  at  (u,v)  and  T  (u,v)  be  the  half-line  of 
R  R 


T  (u,v) \ { 0 }  that  is  tangent  to  W^(u,v),  define  T^(u,v)  similarly. 
R  R  R 


It  is  clear  that  the  evolution  map  of  the  equation  maps  W^(u(R) ,v(R) ) 

R 


to  (u(r) ,v(r) )  i  =  0,1  if  (u(r),v(r))  is  a  solution.  A  similar 


statement  is  true  for  the  action  of  the  tangent  equation  on  Tr. 


Now  let  (u(r)  ,  v(r)  )  be  a  solution  of  (2.5!  i. :  : 

v(r)  <_  0  for  r  <_  R,  some  R,  and  6  (r)  be  the  u-  , ;  :  .*  • 

solution  of  the  angular  variational  equation  (2.  lb').  If 

then  (u (R) ,v (R) )  is  the  upper  half  of  T  (u (R) , v ( R; ; .  b 
R  R 

because  T^fO,!))  points  down  and  if  9(R)  >  C  it  must  : 

00 

the  horizontal  at  least  once  if  (u(R)  ,v(R) )  were  pointir.u  :  -w: 
would  have  passed  through  the  horizontal  twice.  From  the  inn:  : 
in  Theorem  2.2  it  cannot  rotate  through  more  than  an  angle  of  -  - 
but  passing  through  the  horizontal  twice  would  require1  this. 

Now  suppose  (u^(r),v  (r))  and  (u2<r) (r) )  are  two  solar: 
of  the  problem.  Let  u^{0)  <  u2(0)  and  suppose  there  is  no  solutiu 
u(r)  with  u^CO)  <  u(0)  <  u2(0)  nor  any  with  u(0)  <  u ^ ( 0 ) ,  then 
either 

(2.24)  (Ul(R) ,v  (R))  e  W^(u2(R) ,v2(R)) 
or 

(2.25)  (u^R)  .v^R)  )  e  W°(u2(R)  ,v2(R)  ) 

for  all  R  sufficiently  small.  (2.25)  is  impossible  by  comparison 
with  equation  (2.18).  If  (2.24)  holds  then  since  T (u  (R) , v2 (R) }  i 
bounded  away  from  the  horizontal  and  (u  (R) ,v  (R) )  points  up,  fo 
R  sufficiently  small  there  is  a  point  (u,v)  e  W*  (u0(R) ,v9 (R) )  wit 
v  >  0  and  u  >  u  (R)  .  Buth  then  since  (u  (R) , v  (R) )  is  pointin 
down  there  must  exist  another  point  (u^.v^)  with  the  properties 
(1)  v^  <  0  and  (2)  T^fUyV^)  is  the  lower  half  of  T(u^,v^)  whic 
is  a  contradiction  to  the  Remark  on  p.  23. 


Remark:  We  could  have  proved  uniqueness  by  reference  to  the  ful 

PDE  (1.3),  the  spectral  properties  of  Theorem  2.2  mean  that  if  there 
was  more  than  one  solution  their  Morse  indices  would  not  have  added  u 
correctly.  What  we  have  done  is  a  geometric  version  of  the  same  idea 
and  is  done  without  departure  from  the  ODE. 


Chapter  3 


Generalities 

I.  LOCAL  EXISTENCE  AND  UNIQUENESS 

In  this  chapter  we  compile  the  appropriate  generalities  about 
systems  of  reaction  -  diffusion  equations,  everything  is  standard 
material  except  for  section  IV.  Apart  from  some  reversion  to  the 
scalar  case  in  section  II  we  will  consider  the  system  of  equations 
(1.1). 

One  of  the  basic  spaces  (1.1)  can  be  solved  in  is: 

B  =  (u  :  ]Rn  -*■  IR^Iu  is  bounded  and  uniformly  continuous}  , 


where 


To  be  more  precise,  given  u(x)  e  B  there  is  a  unique  solution 
u(x,t)  which  is  a  distributional  solution  of  (1.1)  with  u(x,t)  e  B 
for  sufficiently  small  t  and  u(x,0)  =  u(x). 

The  fundamental  solution  of  the  heat  equation  is 

(3.1)  K(x,  t)  =  (4trt)"n//2  exp[-|xj2/4t] 

It  is  a  standard  matter  to  show  that  u(x,t)  e  C([0,T],B)  (the  con¬ 
tinuous  functions  from  [0,T]  into  B)  is  a  solution  to  (1.1)  with 
u(x,0)  =  u(x)  a^d  only  if  each  component  ui(x,t)  satisfies  the 

integral  equation 

t 

(3.2)  u.(x,t)  =  K(x,t) *u. (x,t)  +  J  K(x,t-s) *f  (u(x,s) )ds 

1  1 


32 


33 


where  *  denotes  convolution  with  respect  to  x  (it  is  here  that  the 
boundedness  of  functions  in  B  enters) . 

If  (3.2)  is  considered  as  a  system  of  equations  (i  =  l,...,m), 
the  right  hand  side  is  a  mapping  on  C([0,T],B).  A  contraction  mapping 
argument  then  shows  that  for  small  enough  T  (3.2)  can  be  solved  for 
any  u(x)  e  B,  and  furthermore  T  only  depends  on  j|u|!  .  An  appli- 
cation  of  the  Gronwall  Inequality  shows  that  the  solution  is  unique  and 
continuous  jointly  in  time  and  initial  data.  For  details  of  these 
arguments  see,  for  instance,  Rauch  and  Smoller  [1]. 

To  sum  this  up,  the  solutions  define  a  local  semiflow  on  B,  i.e. 
if  we  define 

S  (t)  u  (x)  =  u  (x ,  t) 

then  S  is  a  local  semiflow  on  B. 

II.  COMPARISON  PRINCIPES 

To  obtain  qualitative  information  about  solutions  of  a  parabolic 
equation  the  main  tool  is  usually  the  maximum  principle.  We  collect 
here  some  of  the  corollaries  of  the  maximum  principle  that  give  order 
relations  between  the  solutions  of  a  reaction  -  diffusion  equation. 

The  first  three  principles  refer  only  to  the  scalar  equation,  so 
we  assume  u  e  IR  .  Any  functions  u(x,t)  mentioned  are  continuous  on 
some  interval  [0,T]  into  B,  in  the  obvious  fashion.  Since  smooth 
functions  are  dense  in  B,  we  can  assume  u(x,t)  is  smooth  when 
proving  inequalities  such  as  those  below,  as  the  general  case  then 
follows  by  taking  limits. 


G  and 


(3.3) 


■  u 
it 


Au  -  f(u) 


in  some  domain  ft  *  [0,T]  where 


3R  .  Assume 


(3.4) 

(3.5) 


u(x,t) 


u(x,0)  <  0 


then 


(3.6) 


u(x,t)  <0  on 


[0 ,7 J  . 


Principle  3.2.  Suppose 


(3.7) 


-  Au  -  f(u)  >  —  -  Av  -  f ( v ) 
dt  "  ot 


in  some  domain  ft  *  10, T]  where  ft  c  3R  .  Assume  further 


(3.8) 

u(x,t) 

1  v(x,t) 

on 

3ft  x  [ 0  , T] 

(3.9) 

u(x,0) 

1  v(x,0) 

then 

(3.10) 

u(x,t) 

>  V(x,t) 

on 

ft  x  [0  ,  T] 

Proof  of  Principles  3.1  and  3.2.  3.1  clearly  follows  from  3.2  by 

setting  v(x,t)  =  0  and  reversing  the  inequalities.  3.2  follows  fr 
the  standard  linear  maximum  principle  since,  by  the  mean  value  thoo>' 
f(u)  -  f (v)  =  f 1 ( n ) (u-v)  and  so  (3.7)  becomes 

—  (u-v)  -  A  (u-v)  -  f‘  (n)  (u-v)  >_  0 
a  t 


This  is  the  standard  argument,  see  Protter  and  Weinberger  MJ. 


A  function  u(x,t)  that  satisfies  (3.10),  for  any  solutions 


Proof.  Let 

V  (  X  ,  t  )  r  U  (  X  ) 

then 

v(x,t)  is 

a  subsol utior. 

for  (3.12). 

So  u(x,h)  > 

u  (x  ,  0)  ,  if 

u (x , t) 

satisf ies 

(3.12).  Consider  the  rune 

tion  w(x,t) 

=  u(x, t+h) , 

since 

it  satisfies 

wt  = 

Aw  +  f ( w)  , 

by  Principle 

3.2  u  (x , t+h) 

>  u (x 

,  t)  .  u  (x  ,  t) 

is  therefore 

nondecrcas- 

ing  and  if  there  is  a  solution  w(x)  to  (3.13)  we  obviously  must  have 

u(x,t)  <_  w(x)  and  so  lim  u(x,t)  =  v(x)  exists  and  v(x)  <_  w(x)  .  It 

t-*00 

remains  only  to  show  that  v(x'  satisfies  Av  +  f (v)  =  0. 

We  have 

t 

u(x,t+t)  =  K(x,t) *u(x,t)  +  j  K(x , t-s) *f (u (x , s+  i) ) ds 

0 

By  repeated  applications  of  the  dominated  convergence  theorem 

t 

v(x)  =  lim  u(x,t+T)  =  K  (x ,  t)  * v  (x)  +  j  K (x , t-s) *f ( v (x) ) ds 

T-XXD  0 

Since  v(x)  is  then  smooth  it  must  satisfy 

v  =  Av  +  f (v) 

But  v  is  independent  of  t,  so  v(x)  satisfies  Av  +  f (v)  =  0. 

This  proof  is  essentially  that  of  Aronson  and  Weinberger's 
proposition  2.2  [2]. 


We  now  return  to  the  full  system  of  equations  (1.1)  and  define 

invariant  rectangles  for  the  semiflo'w  S(t)  on  B.  Suppose 
m 

p  =  Y  [Vbj  is  a  rectangle  in  3Rm  .  P  is  an  invariant  rectangle 
i=l  1  1 

if  the  set 


M(P)  =  {u(x) 


i  =  1 ,  .  .  .  ,m 


=  (u, (x),...,u  (x) )  a.  <  u.(x)  <  b., 
1  m  1—1  —  l 


and  all  x  <:  )R  } 

is  positively  invariant  (see  definition  1.4).  Weinberger  [1]  and 
Chueh,  Conley  and  Smoller  [1]  have  given  sufficient  conditions  for  a 
rectangle  to  be  invariant,  we  give  these  as  Principle  4. 


Principle  3.4.  A  rectangle  P  = 


II  [a.,b.]  is  invariant  for  the  ' 
i=l  1  1 


semiflow  of  (1.1)  if,  for  all  i  =  l,...,n 


f.  (u)  <0  if  u.  =  b. 

l  —  ii 


f  .  (u)  >  0  if  u  .  =  a . 

l  —  ii 


Remark:  The  condition  says  that  f  does  not  point  out  on  the 

boundary  of  D. 


III.  GLOBAL  EXISTENCE  AND  SMOOTHING 

For  equation  (1.1)  it  was  remarked  in  section  I  that  the  local 
time  of  existence  depends  only  on  !! u  I i  co  anc^  50  an  a  Pri°ri  estimate 
on  ||u(x/t)||  will  suffice  to  prove  global  existence.  If  there  is  an 
invariant  rectangle  P  then  there  is  a  constant  C,  so  that  if 
u(x)  €  M (P)  then  ||u(x,t)||  <  C.  This  means  that  S(t)  induces  a 

global  semiflow  on  an  invariant  rectangle.  From  now,  we  shall  assume 
that  (1.1)  admits  an  appropriate  invariant  rectangle  and  the  semiflow 
will  be  restricted  to  it. 


a. 


Now  suppose  u(x,t)  is  a  smooth  solution  of  i  1 . 1 1  wit: 

u(x,0)  e  M (P)  on  a  set  D  =  *  [0,T]  or  .  -  wr.  r- 

closed.  Let  c  int  D  with  closed,  a  standard  ■  t  m  it 

sup|Vu|  in  terms  of  supjuj  and  information  from  tii*.  •  •  .at  i  . 
D1 

the  author  does  not  know  of  a  proof  in  the  literature,  short  of  • 
Schauder  estimates  which  are  unnecessary  here,  we  sketch  a  ;  r  •  : 
Theorem  3.1.  Under  these  assumptions,  there  is  a  constant  do; 
the  distance  between  and  3D  and  P,  say  C,  so  that 

(3.14)  sup | "u |  <_  C  supjuj 

D1  ° 

where  |lu|  =  £  |?u,|  etc. 
i 

Proof.  Let  £(x,t)  be  a  smooth  cut  off  function  so  that 

if  (x,t)  e  D1  C(x,t)  =  1 
if  (x,t)  e  Dc  £ (x,t)  =  0 

and  £(x,t)  e  [0,1]  always.  Consider  the  system  of  equations 

(3.15)  vk  -  Avk  =  gk(x,t) 


k  k 

on  D,  where  g  is  smooth  and  let  v  be  a  smooth  solution.  Set 


w  =  t  I  lVvJ  2  +  A  I  vi 


Let  K  be  a  constant  (bigger  than  1)  that  satisfies 


K  >  |D0tc|j 


* 

I  am  grateful  to  Professor  L.  C.  Evans  for  teaching  me  this  method, 
is  also  related  to  an  estimate  in  Chueh,  Conley  and  Smoller  [1] . 


where  a  is  a  multi  index  fa.  ^2  and  1  j  2 .  At  a  point  cf 
maximum  interior  to  D  we  must  have 

wfc  -  Aw  >_  0 


computing  and  estimating,  we  get 

0  1  &2(w  ~  Aw)  1  52  I  |Vvk|2(6K2  -  2A)  +  2;2A  £  vkgk 

k  k 

+  2?4  I  Vvk*Vgk  +  8K2^2  l  |Vvk|2  +  8A C  l  uk7uk*V; 
k  k  k 

where  the  last  two  terms  are  arrived  at  by  using  the  fact  that  w 
for  all  i.  For  the  last  term 


?vkVvk-V?  =  (£Vvk)  •  (vk\7£)  < 


;21 Vvkl 2  +  a (vk) 2  I  V C 1 


2a 


for  any  a  >  0,  so  let  a  =  4,  then 

k,2 


0  <1*1  \Vv*\Ul4KZ  -  A)  +  2CZA  l  vV  +  2X?  l  7vK-Vg* 
k  k  k 


+  16AK  ^(vk)2  . 

k 


k  k 

Now  let  g  (x,t)  =  f,  (u(x,t))  and  v  =  u,  ,  since  u(x,t)  t  M 

k  k 


for  all  t  >  0  and  f  is  smooth  it  can  be  seen  that  there  are  con¬ 


stants  L,  and  L„  so  that 
1  2 


l  Vk-li  l  K 


l  Vu  -Vf  (u(x,t))  <  L  l  |Vu 


Plugging  these  into  the  above 
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0  <  y  | Vu  I  2 ( ( 14  +  2L  )K2  -  X))  +  (2KXL  +  16XK) 7  ! u,  '2 
—  ,  k  2  1  "  k 

k  k 


Taking  X  large  enough  the  quantity 

2KLX  +  16XK 
X  -  (14+2L2)K2 


is  positive  and  so 


This  is  true  at  a  point  of  maximum,  so  throughout  D 


5  l  l^uj  +  X  l  |u  |  <  (c+X)  sup  l  ju J' 

k  k  D  k 


but  in  D  C,  =  1  and  (3.14)  follows. 


As  an  application  we  can  set  D  =  Rn  *  [0,°°)  and 
=  3Rn  x  [t0,“)  .  If  u(x,t)  is  a  smooth  solution  of  (1.1)  we  get  a 

bound,  on  sup | Vu (x , t) [  for  all  t  >_  tQ  since  u(x,t)  stays  in  M. 

_n 

]R 

Suppose  we  replace  the  sup-norm  topology  on  B  with  the  compact- 
open  topology,  then  this  says  that  M  •  [t,®)  is  precompact  in  M  for 
each  positive  t,  where  M  consists  of  smooth  functions  in  M.  Since 
clQ (M)  =  M  and 

cl  (M)  •  [t,°°)  c  cl.(M  •  [t,®))  c  cl  (M  •  [t,®)) 
u  u  c 
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where  cl  and  cl^  refer  to  the  sup  norm  and  compact  open  topologies 
respectively,  M  •  [t,°°)  is  precompact  in  M,  with  the  compact-open 
topology . 

To  complete  this  picture  we  must  show  that  the  topology  can  be 
switched  without  destroying  the  semiflow  property.  The  only  property 
to  be  proved  is  continuity,  this  will  be  done  in  the  next  section. 

IV.  CONTINUITY  IN  THE  COMPACT-OPEN  TOPOLOGY 

Given  an  invariant  rectangle  P  for  the  semiflow  S(t)  on  B, 
when  referring  to  the  associated  positively  invariant  set  M(P)  c  b  we 
shall  assume  it  is  endowed  with  the  compact-open  topology. 

In  this  section  we  show  that  if  P  is  an  invariant  rectangle  the 
mapping 

S  :  [0,”)  x  M (P)  -*  M(P) 

given  by  S(t)u(x)  =  u(x,t)  is  continuous.  We  could  not  replace  the 
topology  on  all  of  B  and  get  a  continuous  mapping.  The  a  priori 
bound  implicit  in  M  is  essential. 

Suppose  G  =  {G^}  is  a  sequence  of  compact  subsets  of  ®n  so 
00 

that  I*n  =  bJ  .  then  the  following  quantity  is  a  metric  on  M 

00 

d  (u,v)  =  £  —  sup  |u(x)  -  v(x)  | 

U  i=0  21  xeG. 

l 

To  prove  continuity  of  S  on  [0,°°)  x  m  it  suffices  to  prove 
continuity  of  S(t)  on  M  which  is  uniform  with  respect  to  t  in  some 


compact  interval.  This  follows  from  the  fact  that  for  fixed  ; 
S(t)u  is  continuous  in  t,  from  section  1. 

We  will  find  a  sequence  of  compact  sets  G  =  {G. •  so  that 
0  <  t  <  T  we  have  an  estimate  of  the  form 


(3.16)  d  (u(x,t) ,v(x,t) )  <  c(T,N)d  (u(x) ,v(x) )  +  :(T,N) 

Cj  Cj 


with  N  a  positive  integer  and 

6  (T,N)  -»  0  as  N  -*■  00 

\ 

c  (T  ,N)  -*■  00  as  N  -*•  00 

Such  an  estimate  clearly  performs  the  desired  function  and  we  prove  it 
in  the  following  theorem. 

Theorem  3.2.  In  the  above  notation  S  :  [0,“)  *  M  -*  M  is  continuous. 
Proof.  Recalling  the  heat  kernel  K(x,t)  from  section  I,  with  an 
abuse  of  notation  we  think  of  K(x,t)  as  an  m  *  m  diagonal  matrix 
with  this  kernel  as  each  diagonal  element 


(3.17)  u(x,t)  =  /  K(x-y,t)u(y)dy  +  /  /  K (x-y , t-s) f (u (v , s) ) dyds 
,n 


3R 


0  _n 
]R 


Let  A  c  B  both  be  compact  sets  and  use  the  notation 


|gU) lA  =  sup  |g(x) |  =  sup  £ |gi(x)  | 
xeA  xeA 

|u(x,t)  -  v(x,t)  I  £  |/  +  J  K  (x-y ,  t)  [u  (y)  -  v(y)]dy'| 

n  A 

B  BC 
t 

+  /  |/  +  /  K(x-y,t-s) (f (u(y,s) )  -  f (v(y,s) ) )dy |  ds 
0  B  B°  A 


<  1 u ( x )  -  v(x) |  +  ce(A,B,t)  + 

B 


J  (k|u(x,s)  -  v(x,s)|  +  ce (A,B, t-s) ) ds 

0  B 


where 


and 


c  >  {sup  |u(x)  j  :  u(x)  f  M '■ 


k  >  {sup  D  f(u)  :  1  <  1  <  m  and  u  i  Hi 
1  u .  1  —  — 

i 


e(A,B,t)  =  J  |K{x-v,t)  jftdy 
BC 


Now  let  the  {G.}  be  concentric  balls  of  radius  R..  We  use  the  nota- 
l  i 

oo 

tion  G (n)  for  the  sequence  of  sets  {G.  „  and  d„  ....  for  the 

X+N  1=0  G  (N) 


associated  metric,  then 
(3.18)  d 


G(N) (U,V)  -  2  dG(u-v)  • 


For  each  i,  substitute  A  =  G.  and  B  =  G.  and  sum  over  i  with 

l  l+N 

the  appropriate  weighting. 

d  (u(x,t),  v (x , t) )  <  d  (u,v)  +  c  l  ■—  e(G. ,G.  ) 

G  —  G(N)  .  l  l+N 


t  . 

+  J  (k  d  .  (u(x,s)  ,  v(x,s) )  +  c  l  — 7  c(G. ,G  ) )ds 
‘  G(N)  .  .1  l  l+N 


i  2 


where  e(A,B)  =  e(A,B,T)  and  we  have  used  the  fact  the  e(A,B,t>  is 


increasing  in  t.  Using  (3.18) 
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d  (u(x,t),  v (x  1 1) )  <  2Nd  (u,v)  +  c(l+t)  I  —  e(G.,G  ) 

vj  —  Vj  ^  1  1  +  IM 


N  r 

+  k  2  J  d  (u(x,s) ,  v(x,s))ds 
0  G 


By  Gronwall's  Inequality 


N  N 

w  1-2  T  k2  T  1 

d  (u  (x ,  t) ,  v(x ,  t)  )  <  2  e  d  (u,v)  +  c(l+T)e  E  — r  e (G. , G  ) 
G  —  G  .  _i  l  l+N 

l  2 
N 

N  k2  T 

In  the  notation  of  (3.16)  c(T,N)  =  2  e  and  5(T,N)  = 

N 

k2  T  1 

c(l+T)e  E  — r  e(G.,G.  It  remains  to  show  that  as  N  -*  °° 

.  -l  l  l+N 

l  2 

6(T,N)  ->■  0  for  fixed  T. 


eCG^G.^)  =  (4ttT)  n^2  j  sup  exp(-  |  x-y  |  ^/4T)  dy 


_c  xeG . 
Gi+N 


recalling  that  G^  is  a  ball  of  radius 


=  (4itT)  n /  exp  (-  ( |y  |  -RJ  Z/4T)  dy 


G? 
i+N 


( 4ttT)  /  rm  iexp{-  (r-R  . )  ^/4T)  dr 

R  1 

i+N 


m-1 


=  (4ttT)  n  /  ( q+R . ) "exp(-q“/4T)dq 

R.  -R.  1 

i+N  l 


It  is  a  simple  matter  to  obtain  estimates  of  the  form 

/  qSe  Cq  dq  <_  p(R)e  CR 
R 


5 


If  m  * 


1,  and  if  m  =  1  it  is  even  simpler. 
i+N+1 


Dividin'?  the  •  ' 


in  parentheses  by  3 


we  reach  the  desired  conclusion. 


I.  MAIN  CONTINUATION  ARGUMENT 

The  construction  given  in  section  II  will  allow  us  to  pass  struc¬ 
ture  between  different  invariant  subflows  of  the  semiflow  for  the 
reaction  -  diffusion  equation.  These  subflows  will  be  related  to 
symmetries  of  the  equation,  whence  their  invariance.  In  particular  we 
shall  pass  from  one  space  dimensional  behaviour  to  spherically  symmet¬ 
ric  behaviour. 

We  firstly  express  this  in  the  abstract  setting  of  a  compact  semi¬ 
flow  on  a  fiber  bundle. 

Definition  4.1.  A  semiflow  S(t)  on  a  space  Y  will  be  called  compact 
if  Y  •  [t,°°)  is  precompact  in  Y  for  every  t  >  0. 

We  shall  consider  a  compact  semiflow  on  a  fiber  bundle  but  the 
structure  to  be  continued  will  lie  entirely  within  a  certain  subspace. 
For  the  purpose  of  continuation  the  relevant  features  of  the  subspace 
are  given  by  the  following  definition. 

Definition  4.2.  Suppose  it  :  E  -*  A  is  a  fiber  bundle.  A  subspace 
F  c  e  is  said  to  be  continuous  if  (1)  F  is  closed  and  (2)  tl  is 
an  open  map. 

Note  that  any  fiber  bundle  is  a  continuous  subspace  of  itself.  If 

E^  =  ft  "'■(X),  continuation  involves  deducing  behaviour  in  the  fibers 

E  from  that  in  E.  ,  if  X  is  close  to  X  .  If  A  =  E  n  A,  that 
A  Ag  0  A  X 

A  is  continuous  guarantees  that  this  can  be  done  from  A  to  A  . 

A  A 


47 


■“»  7s 


The  pieces  of  structure  we  shall  continue  are  isolated  invariant  sets 
and  attractors  (see  Conley  [1]  and  Yung  [1]}. 

Let  S(t)  be  a  semiflow  on  Y  and  N  be  any  closed  subset  of  Y, 
define  the  closed  set 

N°=  {y  e  N  :  y  •  [O,®)  c  n} 

A  set  Icy  is  an  isolated  invariant  set  if  I  c  int(N)  for  some 
closed  set  N  and  ui(N^)  =  I,  N  is  then  called  an  isolating  neigh¬ 
bourhood  for  I.  An  isolated  invariant  set  I  is  an  attractor  if  it 
possesses  an  isolating  neighbourhood  N  with  the  property  that 
u(N)  =  I,  such  an  N  is  called  an  attracting  neighbourhood. 

It  is  the  property  of  being  an  isolating  (or  attracting)  neigh¬ 
bourhood  that  is  stable  under  perturbation  not  that  of  being  an  iso¬ 
lated  invariant  set  or  an  attractor.  If  N  is  an  isolating  neighbour¬ 
hood  we  call  the  set  it  isolates  I(N).  Even  though  N  remains  an 
isolating  neighbourhood  for  nearby  semiflows,  I(N)  may  change. 

Now  suppose  we  have  a  compact  semiflow  S(t)  on  the  total  space 
E  of  a  fiber  bundle  and  that  each  fiber  is  invariant.  If  F  c  E  is 
invariant,  that  it  be  a  closed  subspace  of  E  means  the  induced  semi¬ 
flow  on  F  is  also  compact.  In  the  following  theorem,  N  is  a  closed 
set  with  non-empty  interior. 

Theorem  4,1.  Suppose  A  is  a  continuous  invariant  subspace  of  a 

bundle,  which  carries  a  compact  semiflow,  if  N  =  N  n  A  is  an  iso- 

A  A 

0  0 

lating  (attracting)  neighbourhood  in  A  then  N  =  N  n  A  is  a  non- 

A  A  A 

o 

empty  isolating  (attracting)  neighbourhood  in  A^  for  X  sufficiently 
close  to  XQ. 


i 


is  an  isolating  neighbourhood  then  there  exists  a 


Proof .  If 

0 

closed  set  K,  and  an  open  set  V,  both  in  A,  such  that 

X  „  c  X  a  , 


oj(N  )  c  lnt(K  )  c  v  c  n 
A0  A0  A0  A0 

But  then  there  is  a  closed  K  c  E  with  int(K)  j-  £  and  an  open  set 

V  c  E  so  that  K  =  K  n  A  and  V  =  V  n  A  ,  so  that  K  V  ■-  N 
A0  A0  A0  0 


(4.1)  u(N  )  c  Int  (K  n  A  )  c  v  r,  A  c  N 

A0  A0  A0  Ao 

Let  N  =  N  n  A,  since  A  is  continuous  n  tt  ^  ( X )  is  non-empty 

for  A  sufficiently  close  to  A..  Now  N  n  tt  ^(A)  =  N  n  A,  so  it 

0  A  > 


suffices  to  show  that 


w  (N°  n  tt  1  (B)  )  c  Int  (K  n  A) 
A 


for  some  neighbourhood  B  c  A  of  XQ. 

Suppose  (4.2)  were  not  true  then  there  would  be  a  decreasing 

sequence  of  closed  sets  B  so  that  n  B  =*  X-  and  (4.2)  were  false 

n  n  0 

for  B  .  For  each  n,  pick  ay  in  w(N^  n  tt  ^(B  )) 
n  n  An 

n  (N  \Int(K  n  A) )  .  Since  the  semiflow  is  compact  on  A,  w(N°  n  tt  ''"(B,)) 
A  A  1 

is  compact  and  since  A  is  closed  in  E,  w(N^  n  tt  ^(b  )) 

n  J_ 

n (N  \Int(K  n  A) )  is  also  compact.  But  then  the  set  {y  }  has  a  limit 
a  n 

point  y  and  by  the  product  structure  yet  *(A  )  and  clearly  also 

y  €  N  n  Int(K  n  A)  which  implies  y  e  N.  \Int(K  n  A  ) .  But 

A  0  -i  A°  A° 

y  e  w  (N  n  tt  (B) )  n  A.  .  It  is  not  hard  to  see  that  w(NU  n  tt  i(B))  = 
A  Aq  A 

U  (N  n  tt  (A))  and  so  y  e  w(N^  n  tt  1(A  ))  =  w(N^  )  but  this  is  a 
AeB  A  AO  '  0 

contradiction  to  (4.1). 


The  same  argument  proves  the  theorem  for  attract!: 


hoods,  one  only  need  go  through  the  proof  removing  all  t h*. 
of  0 . 


A  much  stronger  statement  is  really  true  here,  namely  :  ■  • 

generalized  Morse  index  (see  Conley  fl],  Yuno  [1!)  should.  ;o:h  mu- 


through  the  fibers  of  such  a  subspace  but  that  the  above  two 
structure  continue  is  sufficient  for  our  purposes. 


II.  GENERAL  APPLICATION 


In  this  section  we  will  construct  a  product  space  (i.e.  a  trivial 
fiber  bundle)  which  portrays  both  the  spherically  symmetric  behaviour 
in  (1.1)  and  its  one-dimensional  behaviour  and  further,  allows  continu¬ 
ation  between  the  two. 

We  start  with  the  product  space  E  =  M  *  10,1]  and  extend  the 
semiflow  on  M  by  leaving  the  second  co-ordinate  fixed.  We  will 
describe  the  continuous  subspace  A,  that  is  of  interest,  by  its 
fibers. 

If  the  underlying  spatial  domain  in  (1.1)  is  ]Rn  then  fix  a  co¬ 
ordinate  system  x  =  (x.  , . . .  ,x  )  on  it.  VJe  describe  A,  =  A  M  ' 

In  \ 

in  two  different  cases. 


X  <  1  u  e  A,  if  u  is  spherically  symmetric  with  the 
A 

origin  considered  as  being  at  (t“"»  0,...,0) 

-L  “  A 


?  =  1  u  r.  A._  if  u  only  depends  on  . 

Note  that  as  X  -*  1,  the  origin  approaches  ~r  alone  the 

x. -axis  and  so  the  functions  in  A,  tend  to  those  in  A,,  this  is 
1  a  1 

what  makes  A  a  continuous  subspace  and  so  is  the  essential  content  of 

Proposition  4.1.  All  the  functions  in  A.  are  functions  of  essentially 

A 

one  independent  variable,  if  \  <  1  this  is  the  radial  variable,  if 
>■  =  1  it  is  x^.  For  each  )  <  1  let 

i 

X  X 

x  =  r  ~  W 

where  is  the  radial  distance  from  the  oriain  at  (r—-,  0,...,0). 

1”  A 

If  X  =  1  let  =  x^.  Considered  as  a  function  of  x,  for  every 
we  have  that  -  x^  on  the  x^-axis.  Given  any  u(x)  ■  A.  it  is 
clear  that  it  can  be  seen  as  a  function  of  x^  so  we  often  write  u(x) 
as  u(x  )  . 

Proposition  4.1.  A,  as  described  above,  is  a  continuous  invariant 
subspace  of  M  x  [0,1]  [0,1]. 

Proof .  A  is  obviously  invariant  since  A  is  invariant  under  rigid 
motions.  We  need  to  show  that  n  :  A  -*•  [0, 1]  is  an  open  map  and 
that  A  is  closed. 

To  show  that  nj^  is  an  open  map,  an  open  set  C  around  u  has 
the  following  lengthy  description 

C  =  {(v,X)  :  |v'x)  -  u(x) |  <  €  for  x  e  K,  K  compact,  v  c 
X  e  I  where  I  c  [0,1]  is  open}. 


it  must  be  shown  that 


C  intersects  all  other  fibers 


If  u  e  A. 

A0 

for  X  sufficiently  close  to  >  This  is  clearly  only  a  problem  if 

0 

Aq  =  1.  Suppose  u  •_  ,  then  u  is  a  function  of  x^  alone,  say 

\  1 

u  =  u(x^) .  Then  consider  v,  •  A,  qiven  by  (x  )  =  u(x  )  if 

XX  X  X 

x  >  -  - — — ,  Each  v  and  u  agree  on  the  half  line  x  >  -  ■ — r  so 

it  is  clear  that  for  a  given  compact  set  K,  X  can  be  chosen  suffi¬ 
ciently  close  to  1  so  that  jv  (x)  -  u(x)  j  <  c  for  x  >:  K. 

To  show  that  A  is  closed,  again  we  need  only  worry  about 

u  (x)  e  A,  with  X  1 .  We  must  show  that  there  is  a  u(x)  c  A, 

n  X  n  1 

n 

such  that  u^(x)  "*■  u(x)  uniformly  on  compact  sets  if  u^(x)  converges 

in  M.  Restricting  each  u  (x)  to  the  x-axis  we  get  a  seauence  of 

n  1 

functions  u^(x^)  that  converge  to  a  function  u^fx^) ,  we  define 

u(x)  =  u^(x..)f  then  it  is  obvious  by  the  nature  of  A  that 
X  A 

n 

u^lx)  ->•  u(x)  in  this  topology.  This  completes  the  proof. 


This  seemingly  innocent  and  simple  result  has  fairly  strong  conse¬ 
quences,  for  it  says  that  attractors  will  continue  from  A  to  A 

1  A 

where  X  <  1.  In  A^  we  have  effectively  one-dimensional  behaviour, 
while  in  A^  for  X  <  1  we  have  spherically  symmetric  behaviour.  So 
if  there  is  an  attractor  for  the  one-dimensional  equation  we  automati¬ 
cally  obtain  a  corresponding  attractor  for  the  spherically  symmetric 
equation.  In  other  words,  some  stable  behaviour  in  the  one-dimensional 
case  yields  associated  stable  behaviour  for  the  spherically  symmetric 
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There  is,  however,  a  drawback  as  we  have  on  A^  the  compact-open 
topology,  which  is  quite  inappropriate  for  discussing  stability.  The 
remedy  for  this  is  actually  restricting  the  discussion  to  another 
invariant  subspace  of  A  in  which  stable  behaviour  does  occur.  As  we 
shall  see  in  the  next  section,  in  the  application  to  the  bistable  equa¬ 


tion  the  topology  restricted  to  the  invariant  subspace  will  be  close 
enough  to  the  sup-norm  topology  to  render  an  attractor. 

So,  if  we  have  an  attractor  for  the  one-dimensional  equation  the 
above  procedure  pulls  one  in  from  "infinity"  for  the  spherically  sym¬ 
metric  equation  but  it  tells  us  nothing  about  the  structure  of  the 
attractor  except  that  it  is  the  maximal  invariant  set  in  some  given 
neighbourhood.  For  a  particular  situation  the  problem  then  is  to 
describe  this  maximal  invariant  set,  for  it  contains  the  asymptotic 
information  about  spherical  propagation.  In  the  next  section  we  shall 
do  this  analysis  for  our  standing  example,  the  bistable  equation. 


III.  SPECIFIC  APPLICATION 


We  consider  again  the  bistable  equation 


ufc  =  Au  +  f (u) 


where  u  e  B  and  f  satisfies  (HI) .  We  base  the  set  M  on  the  tri¬ 
vial  rectangle  P  =  [0,1] ,  which  is  easily  seen  to  be  invariant.  Con¬ 
sequently  (4.3)  generates  a  global  compact  semiflow  on  M  which, 
collecting  the  pieces,  is  given  by 

M  =  {u  :  ]Rn  ->  )R  |  u(x)  is  uniformly  continuous  and  0  £  u(x)  <  1} 


The  construction  of  section  II  of  this  chapter  autcma:  .  -=1- 

supplies  us  with  the  continuous  subspace  A  of  M  <  [0,1]. 

remarked  at  the  end  of  section  II  to  find  an  attractor  in  A,  •■*••• 

1 

restrict  the  space  further.  In  particular  we  require  that  a.  : 

u  €  be  nonincreasinq  in  >;\  which  is  an  invariant  ecr.dit: 

Proposition  4.1  was  not  proved  for  this  A  but  the  proof  r-:  pain  ' 
only  trivial  modifications  to  cover  this  case,  so  A  is  a  continue  :  • 
subspace  of  M  x  [0,1] . 

For  the  one-dimensional  version  of  (4.3) 

(4.4)  u.  =  u  +  f (u) 

t  XX 

we  know  there  is  a  travelling  wave  u(x-ct)  with  u(-*)  =  1,  u(+“)  - 
and  c  >  0,  which  is  monotone  decreasing  and  so  lies  in  A1  .  Becauo-. 
of  the  topology  in  A^  this  solution  appears  as  a  compact  curve  in  A  ^ 
running  from  0  to  1, 

0 

Figure  4.1 

let  us  call  this  compact  set  W. 

In  this  section  we  shall  apply  a  theorem  that  is  proved  in  chapter 
5  which  says  that  W  is  an  attractor  in  A^  and  analyse  its  continua¬ 
tion  to  . 


Theorem  4.3.  With  respect  to  the  space  A^,  u  ^  C  is  an  isolate..: 
invariant  set,  u  5  1  and  the  set  W  are  both  attractors.  Moreover, 
for  any  z  <  a  and  any  a  •:  ]R  the  set 

(4.5)  U (a, t)  =  {u(x1>  £  A1  |  u (a)  £  z) 

•k 

is  an  isolating  neighbourhood.  There  exists  a  s  <  1  so  that  for  an 

★ 

be®  and  1  >  6  >  6 

(4.6)  V  (b ,  6)  =  { u  (x  ^ )  e  |  u(b)  >_  6; 
is  an  attracting  neighbourhood  of  u  =  1.  Finally  if 

W  c  u(a,e)  u  V(b,<5) 

then  U  U  V  is  an  attracting  neighbourhood  of  W. 

(Notation:  as  above,  the  dependencies  of  U  and  V  will  some¬ 

times  be  suppressed.) 

Proof .  This  follows  from  Lemmas  5.1,  5.2  and  Theorem  5.2  in  chapter  5 
section  2. 

Let  us  now  plug  this  information  into  A  and  analyse  the  continu 
sets.  Consider  the  following  three  subsets  of  M  x  [0,1] 

U(a,e)  =  {(u(x),A)  |  u(a,0,...,0)  <_t,  \  e  [0,1]} 

^(b^)  =  {(u(x),A)  |  u(b,0,...,0)  >_  6,  \  e  [0,1]} 

N  =  U  u  V  . 

Then  U  n  A1  =  U(a,e),  V  n  =  V(b,<5)  and  N  n  A^  =  U  u  V .  So  if 
a,  e,  b  and  6  satisfy  the  hypotheses  in  Theorem  4.3  we  can  conclude 
by  Theorem  4.1  that  U  n  is  an  isolating  neighbourhood,  V  n  A^ 

and  N  n  A  are  attracting  neighbourhoods  for  A  sufficiently  close 

A 

to  1.  We  can  clearly  write 


/ 
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A  n  U(a,£)  =  {u(xS  e  A  |  u(a)  £  e} 

A  A 

n  U(b,6)  =  {u{x'S  e  hy  |  u(b)  >_  &  '• 

As  long  as  X  <  1  we  are  already  amongst  spherically  symmetric  func¬ 
tions  and  so  the  above  statements  about  U,  V  and  N  can  be  trans¬ 
lated  into  ones  about  functions  u  :  Pn  -*•  1R  which  depend  only  on  r. 

We  set 

U  (a,e)  =  {u (r)  |  u(a)  <  e} 

s  — 

V  (b,6)  =  (u(r)  I  u(b)  >  6} 
s  - 

if  a  and  b  are  large  enough  Us  is  an  isolating  neighbourhood  and 

V  is  an  attracting  neighbourhood,  as  is  U  u  V  if  b  -  a  is 
s  s  s 

sufficiently  large.  We  assume  that  these  conditions  on  a  and  b  are 
satisfied  and  analyse  the  invariant  sets. 

Proposition  4.2.  I(Vs(b,6))  =1  if  6  is  close  enough  to  1  and  b 

is  sufficiently  large. 

Proof.  One  can  check  easily  that  under  these  circumstances  one  of  the 
comparison  functions  in  Aronson  and  Weinberger  [2],  see  our  theorem  1.1, 
will  be  less  than  every  element  of  Vg(b,6). 

Since  the  semiflow  is  compact  and  Ug  u  Vg  is  connected, 

o)(U  U  V  )  is  connected  and  since  0  e  U  ,  I(U  )  d> .  The  usual  argu- 
s  s  s  s 

ment  then  shows  that  co(U  u  V  )  contains  I  (U  ),  1  and  orbits 

s  s  s 

running  from  the  former  to  the  latter.  In  particular  I(U  )  is  not  an 

s 


attractor,  however  we  do  have  the  following  proposition;  recall  that 
Aq  is  just  the  set  of  spherically  symmetric  functions  which  are  non¬ 
increasing  in  r. 

Proposition  4.3.  u  5  0  is  an  attractor  in  A  . 

Proof .  It  is  a  standard  maximum  principle  argument  that  u  =  0  is  an 
attractor  in  the  sup-norm  topology.  Take  any  compact  set  K  c  ]Rn 
which  contains  the  origin,  if  functions  in  A^  satisfy  an  estimate  on 
K  they  satisfy  it  everywhere  and  so  a  neighbourhood  in  the  compact 
open  topology  is  contained  in  a  sup-norm  neighbourhood  which  0 
attracts . 

This  proposition  shows  that  I(U  )  contains  something  other  than 

s 

just  u  =  0.  However  we  can  show  that  the  flow  on  IdM  induced  by 
the  semiflow  of  the  equation  is  a  gradient  flow,  that  is  there  is  a 
functional  V(u)  which  is  decreasing  along  orbits  except  on  constant 
solutions.  The  functional  V(u)  is  the  standard  energy  functional 

1  ?  u 

(4.7)  V(u)  =  /  (-  I Vu |  -  /  f (s) ds)dx  . 

'  n  2  '  1  ' 

]R  0 

It  is  a  standard  computation  to  show  that  along  orbits 

(4.8)  V(u(x,t) )  =  -  j  u^dx 

dt  J  n  t 

IR 

The  only  reason  (4.7)  is  not  used  more  often  is  that  the  integrals  are 

often  unbounded.  However,  here  we  can  show  that  I(u  )  c  h^K11) 

s 


which  is  in  the  domain  of  V(u)  and  so  it  i 


1 

i 

I 

i 

Theorem  4.4.  I(U  )  c  H^(Rn)  . 

3  1 
Proof .  We  actually  show  much  more,  namely  if  we  define  t:.’  t  j 

U^'={utU  I  u  •  t  e  U  for  all  t  T  I 

s  1  s 

T  T  ^ 

then  U  •  T  is  contained  in  some  set  Z  which  is  closed  and  , 

consists  of  functions  that  are  exponentially  decreasing.  Since 
T 

I(U  )  c  u  •  T  for  all  T  >  0,  elements  of  I(U  )  are  exponential  lv 
s  —  s 

decreasing  and  we  can  then  estimate  the  derivatives  to  prove  the  pro; - 
sition. 

We  will  define  a  supersolution  to  the  equation  u(r,t)  with  the 
following  properties 

(4.9)  u(r ,0)  >_  e 

(4.10)  lim  u(r,t)  <_  C  e  mr  some  constants  C,  m  >  0 
t~x» 

(4.11)  u(a,t)  >_  e  for  all  t  >  0 

Suppose  that  these  three  properties  can  be  satisfied,  then  the  function 

min{o,  u(r,t)}  is  also  a  supersolution  on  the  unbounded  annulus  r  ^  a. 

T 

If  u  e  U  then  u(a,t)  <_  z  for  t  <_  T  and  since  u(r,0)  =  u(r)  £  ; 
for  all  r 

u(r,T)  <_  min{a,  u(r,T)  } 

T  I 

So  define  Z  as  the  set  of  v(r)  e  that  satisfy 


v(r)  <  min{a,  u(r,T)} 


T 

From  (4.10)  ^  2  consists  of  exponentially  decreasing  functions. 

T>0 

To  define  u(r,t)  we  use  a  method  of  Fife  and  McLeod's  [1]  from 

their  Lemma  4.1.  Let  v(x)  be  anv  solution  of  v  +  f(v)  that  is 

xx 

nonconstant  but  lim  v(x)  =  0,  i.e.  the  body  of  the  fish  in  Figure  3, 
x-*-±°° 

that  also  satisfies  v(a)  >  e.  Set 

u(r,t)  =  v(r-p(t) )  +  cekt 

where  c,  k  and  p(t)  are  to  be  chosen.  If  p(t)  >  0  for  t  >  0 
and  p(t)  is  bounded  as  t  -*■  °>,  c  r  and  k  <  0  then  4.9  -  11  are 
satisfied. 

We  only  need  show  that  u(r,t)  is  a  supersolution  if  u  <_  a  and 

only  on  the  half-line  [c,<»)  where  v(c)  =  a  and  v'(c)  <  0.  We 

compute 

Lu  =  at  -  ;rr  -  f(u)  >  ut  -  urr  -  f(u> 

(4.12) 

=  -p ' v '  +  kcekt  -  v"  -  f (v  +  cekt) 

Consider  two  sets  separately  u  c  [0,6]  and  u  e  (6, a]  where,  for 
constants  and 

f' (u)  <  q^  <  0  for  u  e  [0,6] 

v'(r)  <  q^  <  0  for  u  e  (6, a] 

From  (4.12) 

Lu  >_  -p'v'  +  kcekt  -  (f(v  +  cekt)  -  f(v)) 


If  u  £  [0,6] 


If  p 1 (t)  >0  then 


-  kt , 

Lu  >  ce  (k  -  q^) 

so  pick  k  <  0  but  k  >  q^  then  Lu  0.  If  ue  [6, a]  there  is  a 
K  so  that 

_  1c  t 

Lu  >_  -p'q  +  ce  (k  -  K) 

so  we  can  set 

,  c(k-k)  kt 

p.  =  — -  e 

42 

since  k  <  K  and  <  0 ,  p'  ^  0,  we  can  then  set 


p(t) 


ext  +  z 
kq2 


to  get  Lu  >_  0  for  ail  u  e  [0,a]  .  Pick  z  so  that  p(t)  >_  0  for 
t  ^  0,  clearly  p(t)  ->  z  as  t  -*■  ®  and  choose  c  =  t. 

From  this  it  follows  that  I(Ug)  consists  of  exponentially 
decreasing  functions  and  therefore  these  functions  are  in  L^(®n)  .  To 
estimate  |vu(x)|  for  u(x)  e  I(Ug),  take  ft  =  {r  |  r  RQ}  and 

ft'  ={r  [  r>_R^}  with  >  R^  then  from  Theorem  3.1  with 

D  =  ft  x  [0,®)  and  D'  =  ft'  x  [0,®) 

sup  | V  u |  £  c  sup  | u I 
D'  D 

and  so 
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sup  |  V  u  |  <_  c  sup  |  u  ] 

?.  C. ' 

for  any  u  t  I(U  ) .  Since  c  depends  only  on  R_  -  R,  and  iu  is 
s  0  1 

exponentially  decreasing  it  follows  that  | V u |  is  also,  which  proves 
the  theorem. 


That  I(Us>  inherits  the  semiflow  as  a  gradient  flow  means  that 
it  consists  of  equilibrium  solutions  and  orbits  connecting  them.  All 
the  equilibrium  solutions  must  be  ones  that  would  satisfy  the  boundary 
conditions  in  Theorem  3.1. 

If  we  suppose  that  f  satisfies  (H2) ,  then  there  is  a  unique  non¬ 
constant  equilibrium  solution  satisfying  u(+“)  =  0  (Theorem  3.3).  In 
this  case  KlM  consists  of  u  ;  0,  this  equilibrium  solution,  call 
it  u,  and  orbits  connecting  them. 

In  w(U  u  V  )  there  are  orbits  running  from  I(U  )  to  1.  Let 

s  s  s 

u(r,t)  be  such  an  orbit,  then  lim  u(r,t)  c  i (u  ).  if  t  is  a  large 

t-*~ S 

enough  negative  number  we  can  show  that  u(r,t)  is  exponentially 

decreasing  in  r  and  u(r,t)  e  .  Pick  T,  a  large  enough 

negative  number  so  that  u(r,t)  €  U  for  all  t  <  T,  then  there  is  a 

s  — 

sequence  u  so  that  u  •  t  ->  u(r,T)  and  u  •  t  e  U  for  all 

n  n  n  ns 

n  T 

t<t  so  u  •  t  e  Z  but  then  u(r,T)  e  Z  and  the  above 
—  n  n  n  „  „ 

T>0 

follows  as  before. 

A  consequence  of  the  fact  that  u(r,t)  e  if  t  is  large 

negative  is  that  lim  u(r,t)  is  an  equilibrium  point  in  I(U  )  and 

s 

2 

now  this  limit  can  be  considered  in  the  L  norm  or  the  sup-norm. 


—  j'*'  *>■■**' 


Suppose  again  that  f  satisfies  (H2)  then  the  only  :  - 
solution  u(r,t)  can  approach  in  backward  tine  is  u.  Morc.'Vi  :  , 
Theorem  2.3  we  know  the  spectrum  of  the  linearized  problem  arour.  : 
and  so  we  can  apply  the  invariant  manifold  theorems  lr.  infinite  di: 
sions.  For  the  equation 

(4.13)  u^  -  Au  +  f(u) 

the  linearized  operator  around  u  acting  on  v  is 

(4.14)  Av  +  f*(u)v  =  Lv 

We  consider  o(L)  with  respect  to  the  space  L^(K+)  with  the 

measure  rn  1dr,  then  it  is  standard  that  a(L)  03*  and  th.e 

essential  spectrum  is  contained  in  (-<*>,b]  where  b  =  f  *  ( 0)  and  so 

b  <  0  (by  HI).  But  then  from  Theorem  2.3,  o(L)  n  [0,°°)  consists 

only  one  point  X  and  X  >  0.  Then  from  Theorem  5.2.1  in  Henry  [1! 

there  is  locally  at  u  a  one-dimensional  unstable  manifold  which  :o 

sists  of  all  solutions  that  approach  u  in  negative  time. 

From  this  we  see  that  there  are  only  two  solutions  leaving  u, 

one  must  go  to  0  and  the  other  to  1.  So  if  f  satisfies  (H2)  we 

have  a  complete  description  of  the  spherical  attractor  W  ,  oictori 

s 

ally. 


1 


Figure  4.2 


Chapter  5 


Spherically  Symmetric  Solutions  in  a  Moving  Frame 
I.  EXTENSION  OF  CONTINUATION  ARGUMENT 

The  space  A  c  M  x  [0,1]  of  chapter  4,  section  II  inherited  its 
semiflow  from  M  x  [0.1]  which  comes  from  the  equation  in  question. 

In  this  chapter  we  endow  A  with  a  family  of  semiflows  that  do  not 
leave  the  fibers  invariant.  Studvinq  the  behaviour  in  A^  and  its 
relationship  to  the  rest  of  these  semiflows  will  give  much  finer  infor¬ 
mation  about  the  propagating  solutions.  We  will  not  present  this 
application  until  section  III  of  this  chapter,  but  the  construction  of 
the  semiflow  should  serve  as  motivation  for  the  generality  of  the 
theorems  and  propositions  in  section  II.  The  statements  in  this  section 
are  to  be  understood  in  the  context  of  a  general  system  of  reaction  - 
diffusion  equations,  but  the  remainder  of  the  chapter  is  exclusively 
about  the  bistable  equation  (1.3)  where  f  satisfies  (HI). 

Let  us  call  the  basic  semiflow  on  M,  S(t).  We  shall  define  a 
semiflow,  called  H(t),  on  the  space  A  x  [0,°°).  Firstly  we  define  a 
translation  map  on  A, 

T  :  1R  x  a  -*  A 


by  the  formula  (recall  A  c  m  x  [0,1]) 

(5.1)  T (a ,  u(x,,...,x  ) ,X)  =  (u (x, +a,  x  , . . . ,x  ) , A  *  > 

in  i  4  n 


where 

(5.2) 


a  +  (X/,,) 

\  I  _  _ _ 

1  +  a  +  (X/1_x) 
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It  is  clear  that  if  u(x  .  ...,x  )  c  A.  then  u(x,+a,  x^,...,x  )  •  A.,, 
where  is  given  by  (5.2),  so  T  is  well  defined.  Since  transla¬ 

tion  is  continuous  on  M  and  X'  is  a  continuous  function  of  '•  and 
a.  from  inspection  of  (5.2),  T  :  3R  x  a  *  A  is  continuous.  Intu¬ 
itively  T  translates  functions  in  the  x^  direction  and  moves  them 
over  to  the  appropriate  fiber.  Notice,  from  (5.2)  that  A.  is 
invariant  under  T,  as  it  should  be. 

For  c  e  [0,co)  ,  we  use  T  to  define  a  semiflow  H  (t)  on  A. 

c 

Let  (u(x),A)  e  A,  define 

(5.3)  H  (t)(u(x),X)  =  T(ct) (S(t)u(x) ,X)  . 

c 

If  c  >  0,  as  t  -»<*>,  H  (t)  pushes  all  of  A  onto  A  ,  intuitively 

c  1 

it  pushes  the  origin  out  to  -®.  Consequently  the  asymptotic  analysis 
for  this  semiflow  will  depend  on  its  behaviour  in  A^.  Hc(t)  gives 
the  evolution  of  the  equation  but  viewed  while  moving  out  in  a  radial 
direction  with  speed  c.  Restricted  to  A^  ,  the  semiflow  is  that  of 
the  one-dimensional  equation  in  a  moving  co-ordinate  frame,  so  we  will 
see  how  the  asymptotic  behaviour  of  the  spherical  solutions  is  deter¬ 
mined  by  the  one-dimensional  equation  at  different  speeds. 

For  each  t  >_  0,  define  H(t)  on  A  *  [0,°°)  by 

H(t)  ( (u,X)  ,c)  =  (H  (t)  (u,X)  ,c) 
c 

It  is  clear  that  H  defines  a  semiflow  on  A  *  [0,“)  that  contains 
( t)  by  restricting  to  A  x  { c} . 

We  must  turn  now  to  analysing  the  asymptotic  behaviour  of  the  one¬ 


dimensional  bistable  equation  in  various  moving  co-ordinate  frames. 
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II.  ATTRACTOR  IN  ONE  DIMENSIONAL  CASE 

For  the  standing  example  we  assume  again  that  f  satisfies  (Hi) 

(see  chapter  1).  In  A^  there  are  functions  of  a  single  variable, 

call  it  £,  which  are  nonincreasing  in  £  and  stay  between  0  and  1. 

The  semiflow  H  (t)  when  restricted  to  A  is  that  of  the  equation 

c  1 

(5.4)  ufc  =  u  +cu^  +  f  (u) 

where  £  =  x^-ct,  so  this  agrees  with  previous  notation  when  c  =  0. 

The  set  W  is  the  union  of  0,  1  and  all  translates  of  the 
travelling  wave,  which  is  invariant  under  each  of  the  semiflows  H^(t). 

The  main  goal  of  this  section  is  to  prove  that  W  is  an  attractor  for 

*  *  _ 

c  £  [0,c  ],  some  c  >  c,  relative  to  A^ . 

The  strongest  stability  results  for  the  travelling  wave  are  due  to 

Fife  and  McLeod  [1] .  Specialising  some  of  their  results  to  the  present 
case  reads  as  follows. 

Theorem  5.1.  (Fife  and  McLeod) .  If  f  satisfies  (HI)  and 
0  <  u(x)  <^1  is  uniformly  continuous,  x  £  IR  ,  and  satisfies 

lim  inf  u(x)  >  a  and  lim  sup  u(x)  <  a 
x-*-00  x-*-+°° 

then  there  is  a  constant  xQ  such  that 

lim  ||u(x,t)  -  u^(x  -  ct  +  xQ)  11^=  0 

t^-yoo 

where  u^(£)  is  the  one-dimensional  travelling  wave.  Further,  the 
limit  is  uniform  over  a  sup-norm  neighbourhood  of  u^. 


This  could  be  described  as  uniform  asymptotic  stability  . : 
travelling  wave.  Fife  and  McLeod  actually  prove  a  much  stror.  :■  r 
namely  exponential  stability,  but  we  shall  not  use  this. 

Before  proving  that  W  is  an  attractor,  the  local  dynamical  • : 
perties  of  the  constant  solutions  u  =  0  and  u  ;  1  must  be  stu.i: 

In  the  topology  of  A^  a  closed  neighbourhood  of  i.e.  on- 

with  non-empty  interior  containing  0,  is  a  set  which  detenus  on  a 
compact  set  K  and  a  positive  function  g(£)  defined  on  K.  Such  a 
set  is  then  defined  as  (recalling  A^  c  m) 

U(K,g)  -  {u  €  A  I  u(£)  £  g(£)  for  £  e  K; 

A  neighbourhood  of  1  is  given  similarly 

V (K,g)  =  {v  e  A1  |  v(£)  >_  g  ( £)  for  £  e  K} 

If  K  =  {a},  then  we  write  U(a,e)  and  V(a,e)  where  e  =  g(a), 

this  is  consistent  with  the  notation  of  chapter  4. 

Throughout  this  chapter  we  shall  apply  the  comparison  principles 

* 

of  chapter  3,  section  II  to  the  equation  (5.4),  it  is  trivial  to  exte 
each  one  to  cover  this  case . 

Lemma  5.1.  If  e  <  a,  for  any  a,  in  the  c-semiflow  on  A  , 
u  =  0  is  an 

(A)  isolated  invariant  set  if  c  >  c  >  0 

(B)  attractor  if  c  >  c 

and  U(a,e)  is  an  isolating  (in  case  (A))  and  attracting  (in  case  ( 


neighbourhood . 


Proof.  To  deal  with  Case  (A)  define 

0 

U  (a,t)  =  {u  e  U(a,c)  :  u*  t0,“)  c  u(a,e)> 
c  c 

where  " •  ”  refers  to  the  action  of  the  c-semiflow.  From  chapter  3, 

it  suffices  to  show  that  w  (U^(a,e))  =  {0},  where  oj  refers  to  th 

c  c  c 

u-limit  set  in  the  c-semiflow.  This  will  follow  if  we  can  show  that 

given  any  neighbourhood  U  of  0  there  is  a  T  so  that 
0 

Uc(a,e)*c  t  c  u  for  all  t  >.  T.  By  Comparison  Principle  II  it  is 

obviously  adequate  to  find  a  solution  to  (5.4)  so  that  if 

v(£)  e  U°(a,e)  then 
c 

v(£)  £  u(£,0) 

and  u(£,t)  •*  0  in  A^. 

From  theorem  5.1,  if  c  <  c,  then  v(£)  £  U°(a,e)  implies  that 

c 

lim  v(£)  <_  a.  So  we  need  only  construct  a  solution  with  u(?,0)  >_  .i 
if  £  <  a  and  u(£,0)  >_  z  if  £  >_  a. 

Transforming  (5.4)  into  a  system  yields 

w'  =  z 

(5.5) 

z'  =  -cz  -  f(w)  , 


If  c  =  0,  the  phase  portrait  is  in  Figure  2.1,  for  c  >  0  it  is 


Figure  5.1 


The  important  fact  here  is  that  for  any  0  <  e  <  a  the  solution  with 
w(0)  =  e  and  z{0)  =0  must  cross  the  line  w  =  c;  in  the  lower  half 
plane,  let  d  be  the  point  that  gives  w^d)  =  a  and  ztc)  <  0.  Now 
define  u(£)  by 


' 

f  ° 

if 

£ 

< 

a  +  d 

u(0  =  < 

w(£-a) 

1 

if 

a 

+ 

d  <  ^  < 

l  « 

if 

£ 

> 

a 

Then  u(£)  satisfies  the  hypotheses  of  Comparison  Principle  III  and  so 
clearly  u(C,t)  -*■  0  as  t  -*■  00  uniformly  on  compact  sets.  3at  also  if 
5  <  a,  u<£,0)  =  u(£)  =  a  and  u(?,0)  >_  c  if  5  >  a.  So  u(£,t)  per¬ 
forms  the  desired  function. 

Now  suppose  c  >  c,  for  case  (B) ,  to  show  0  is  an  attractor  we 

need  u  (U)  *■  { 0 } .  As  before  we  construct  u(£)  so  that  v(£)  £  U 
c 

implies  v(£)  <  u(£)  #  and  u(£«t)  ->0  in  A  .  But  by  theorem  5.1 
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any  function  in  A^  which  satisfies  (i)  t  <  1  im  u(£)  <  a  and  (ii) 

x-^+ct3 

u(£)  =  1  if  £  a  (of  which  there  are  obviously  many)  will  have 
u(£,t)  -*-0  in  A^,  since  such  a  function  will  obviously  majorise  U, 
the  Lemma  is  proved. 

•k  _ 

Lemma  5.2.  There  is  a  c  >  c  so  that  u  1 ,  in  the  c-semif low- 
on  A^ ,  is  an 

(A)  attractor  if  0  <  c  <  c 

*  — 

(B)  isolated  invariant  set  if  c  >  c  >  c 

and  if  1  >  6  >  a  then,  for  any  b,  V(b,6)  is  an  isolating  neighbour¬ 
hood  in  case  (A).  For  each  c  in  (B) ,  there  is  a  6  so  that  V(b,5) 
is  an  attracting  neighbourhood,  for  any  b. 

Proof.  The  proof  is  very  similar  to  that  of  Lemma  5.1.  In  case  (A)  a 

function  u(£)  below  everything  in  V(b,6)  can  easily  be  found  if 

6  >  a  and  by  Theorem  5.1  u(£,t)  -»■  1  in  A^. 

For  case  (B)  ,  by  Theorem  5.1,  V°  consists  of  functions  v(£) 

c 

* 

which  satisfy  v(£)  >_  a  for  all  £.  If  c  is  not  too  much  bigger 

than  c  the  unstable  manifold  of  (1,0)  intersects  the  line  w  =  a 

_  * 

in  the  lower  half-plane  for  the  system  (5.5)  with  c  <  c  <  c  . 


”1 


Remark:  The  conditions  under  which  Lemma  5.2  gives  an  isolating 

neighbourhood  of  u  =  1  are  a  little  cumbersome,  but  the  following 
three  simple  points  will  help. 

(1)  If  V(b,6)  is  an  isolating  neighbourhood  for  a  given  c  >  c 
then  it  is  for  all  c‘  which  satisfy  c  <  c‘  <  c. 

(2)  If  V(b,6)  is  an  isolating  neighbourhood  for  c  >  c,  then 
so  is  V(b,6')  if  6'  satisfies  1  >  6'  >  6. 

(3)  If  the  solution  to  (5.5)  with  z(0)  =  0,  w(0)  =  5  crosses 

the  line  u  =  a  in  the  lower  half-plane  then  V(b ,6)  is  an  isolating 

neighbourhood  for  the  c-semiflow. 

Before  proving  that  W  is  an  attractor  we  must  consider  some 

* 

slightly  more  exotic  neighbourhoods  of  0  and  1.  For  each  c  e  [0,c  ), 

* 

where  c  is  the  same  as  in  Lemma  5.2,  and  each  e  <  a  sufficiently 
close  to  a  the  solution  of  (5.5)  with  w(0)  =  e  and  z(0)  =  0  must 

cross  the  w-axis  between  a  and  1  in  backward  time.  Let  d  <  0  be 

the  largest  number  such  that  z(d)  =  0  and  set  6  =  w(d).  Let 
K  =  [a,a-d]  and  g  (£)  =  w(£  -  (a-d) )  with  domain  K.  Suppose 

cl  3. 

e-d  <  a,  then  we  will  use  the  sets  U  (K  ,g  )  and  V  (K  ,g  )  given 

c  a  a  c  e  e 

by 

U  (K  ,g  )  =  {u  e  A  |  u(£)  <  w(£  -  (a-d))  for  a  <_  £  £  a-d} 

C  3  3  X 

V  (K  ,g  )  =  (u  e  A,  I  u(£)  >  w(£  -  (e-d))  for  e  <  £  <  e-d} 
c  e  e  1  —  — 


The  picture  is 


e 


e-d 


a 


a-d 


Figure  5.3 

U  (K  ,g  )  consists  of  functions  below  the  right  hand  curve  and 
c  a  a 

V  (K  ,g  )  of  those  above  the  left  hand  curve, 
c  e  e 

* 

For  each  c  e  [C,c  ]  these  sets  are  isolating  neighbourhoods  of 

0  and  1  respectively  as  they  are  contained  in  sets  of  the  form  found 

in  Lemmas  5.1  and  5.2.  Namely,  U  (K  ,g  )  c  U(a-d,e)  and 

c  a  a 

V  (K  , g  )  c  v(e,S) ,  the  fact  that  w  exists  implies  that  5  is  close 
c  e  e 

enough  to  1  for  V(e,6)  to  be  an  isolating  neighbourhood  by  remark 
(3)  after  Lemma  5.2.  The  crucial  property  of  these  sets  is  that  they 
give  us  a  positively  invariant  neighbourhood  of  W. 

♦ 

Lemma  5.3.  For  any  0  <  s  <  o  and  any  c  e  [0,c  ] ,  if  e-d  <  a 
then  the  set 


U  (K  ,g  ) 
c  a  a 


u  V  (K  ,g  ) 
c  e  e 


is  positively  invariant  in  the  c-semiflow. 

Proof.  Let  u(U  c  U  (K  ,g  )  then  we  must  show  that 
-  c  a  a 

u(£,t)  c  U  (K  ,g  )  u  V  (K  ,g  )  for  all  positive  t.  So  there  are  two 
c  a  a  cee 


things  to  prove: 
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(1)  If  u(£,i)  e  U  (K  ,g  )  but  u(£,t)  i  U  (K  ,g  )  for  t  >  t 
c  a  a  c  a  a 

but  close  to  x  then  u(£,t)  e  V  (K  ,g  )  and  (2)  the  same  as  (1)  with 

c  e  e 

the  roles  of  and  reversed.  The  proofs  of  (1)  and  (2)  are, 

not  surprisingly,  almost  identical,  so  we  only  prove  (1). 

Now  suppose  u(£,t)  is  a  solution  of  (5.3)  with  u(£,t)  t-  u  (K  ,g  )  . 

C  3  3 

Since  g  on  K  is  a  solution  of  (5.3)  itself  then  u(£,t)  t  U  (K  ,g  ) 

3  3  C  3  3 

for  all  t  <_  t  unless,  by  comparison  principle  I,  u(£,t)  c  U  (K  ,g  ) 

C  3  3 

for  T  t  <  t  and  u(£,t  )  =  g  (£)  for  £  e  3K  ,  i.e.  £  =  a  or 

0  U  3  3 

£  =  a-d.  Since  u  e  A  is  nonincreasing,  if  u(a,T  )  =  g  (a)  =6  we 

A.  0  3 

have  u (£ , t  )  c  V  (K  ,g  )  . 

0  c  e  e 

Suppose  u(a-d,Tg)  =  w(0)  =  e  then  the  following  inequality  holds 

u^  (a-d,TQ)  +  cu^(a-d,TQ)  +  f(u(a-d,TQ))  £w^(0)  +  cw^(O)  +  f(w(0)) 

the  last  terms  on  each  side  are  equal.  Since  w  (0)  =  0  and  u_  £  0, 
if  the  inequality  were  not  true  we  would  have 

U££(a-d'V  "  W££(0) 

but  this  is  impossible  since  u(a-d,TQ)  =  w(0)  ,  u^  <_  0  and 

u ( £ , Tq )  _<  w(£)  for  £  <  0  and  close  to  it. 

The  above  inequality  implies  that  ut(a-d,iQ)  <_  0  and  so  this 

does  not  provide  an  escape  from  U  (K  ,g  )  .  This  completes  the  proof 

C  3  3 


of  Lemma  5.3. 


We  will  show  that  the  set  U  c  V  of  Lemma  5.3  is  ir.  attra  •-  : 

c  c 

neighbourhood  for  W  for  each  c-semif low,  the  travelling  wav- 

and  1  attached).  But  for  the  sake  of  hygiene,  we  give  a  no;  ••  ••  >. 

described  set  that  is  contained  in  U  L  V  .  In  the  definition  v 

c  c 

and  V  ,  6  =  w(d),  so  <5  depends  on  c  and  setting  ‘  =  ‘(c), 
c 

define 

6  =  sup{6(c)  :  c  £  [0,c  ] : 

which  is  clearly  less  than  1.  With  b  =  e-d,  if  &  >_  6 
U (a, e)  u  V  (b ,  6)  c  u  (K  ,g  )  t  V  (K  ,g  ) 

C  cl  ci  C  c:  e 

* 

for  all  c  which  satisfy  0  <  c  <  c  .  Note  also  that  if  b  is 

sufficiently  smaller  than  a  then  W  c  U(a,£)  U  V(b,6).  We  are  free 

to  choose  e  <  a  as  this  is  free  to  be  chosen  in  Uc  and  Vc,  but  : 
must  satisfy  6  >_  6,  which  depends  on  e.  We  will  see  in  chapter  6 
that  it  is  useful  to  have  no  restriction  on  E,  the  restriction  on 
is  unimportant. 

* 

Theorem  5.2.  If  c  e  [0,c  ]  then  W  is  an  attractor  in  the  c-semif low 

and  there  exists  a  6  >  a  so  that  for  any  e  <  a  if 

W  c  u (a , e)  U  V(b, 5) 

then  U  u  V  is  an  attracting  neighbourhood  for  W  in  each  of  these 
semi flows. 

Proof.  Consider  u  (U  U  V)  =  D  since  U  u  V  c  u  u  V  (the  depend- 
-  c  c  c  c 

encies  on  K  ,a  etc.  are  being  suppressed,  it  is  assumed  they  satisfy 

the  requirements  in  the  preamble  for  the  theorem)  we  have 

DC(!)(UuV)cuuV,  by  Lemma  5.3.  0  e  D  and  1  t  D  . 

cccccc  c  c 


Suppose  a  solution 

u(£,t)  is 

in  D  , 
c 

then 

it  is  defined 

for  all 

t  €  K  and  u(£,t) 

tr\ 

a 

o 

c 

< 

o 

for  all 

t. 

(Recall  that 

=  xrc'--) 

Consider  the  case  c  <  c. 

either 

(A) 

u(C,t)  -*  0  as 

t  •>  or 

(B)  u(£,t)  -*•  1  as  t  If  (A)  happens  we  cannot  have 

lim  u(£,»t)  =  0  unless  there  exists  a  t  such  that  u(^»t)  U  ,  l .  <:  ■ 
t-*— 00 

u(£,x)  e  V  but  then  lim  u(£,t)  =  1  since  V  is  an  attractina 
c  c 

p-XJO 

neighbourhood  of  1,  which  contradicts  (A).  So  (A)  implies  that 
lim  u(C/t)  =  1,  which  is  impossible,  so  (A)  is  impossible  and  (B) 

t->— co 

holds.  We  cannot  have  lim  u(£,t)  =1  as  it  is  an  attractor,  so 

£->-—co 

lim  u(£,t)  =  0. 

£—>—oo 

By  analogous  arguments,  if  c  >  c,  we  must  have  lim  u(£,t)  =  0 

t~H-oo 

and  lim  u(£,t)  =  l. 

£->—00 

Since  the  same  statements  obviously  hold  for  wc(uc  u  '  a  con~ 
sequence  is  that  for  c  ^  c 

10  (U  U  V)  =  to  (U  U  V  ) 
c  c  c  c 

The  inclusion  from  left  to  right  is  trivial  but  if  u(£,t)  e  to  (u  IV) 

c  c  c 

for  all  t  then  lim  u(£,t)  =  0  (if  c  <  c)  and  since  U  is  a 

CO 

neighbourhood  of  0,  u(£,t)  must  be  in  w  (U  u  V)  as  well,  it  is  a 

c 

similar  argument  for  c  >  c. 

Suppose  v(£)  e  D  then  clearly  v(£+a)  will  also  be  in  D 

C1  C1 
lim  v(5,t)  =  lim  v(£+a,t) .  Also  there  is  a  v  (£)  e  D  c  u  tv 

t— 00  **—  T  C1  C1  C! 

such  that 


rjLMU  V\***i*Wv 


and 


but  then  ( £  -  (c^-c^JT)  < 

H  ( T ) V  ( £  -  (c  -c  ) T)  =  T( (c  -c, )T)H  (T)V(£  -  (c--c.)T) 
c  T  21  21c,  T  21 

2  1 

=  Hc  (T)T(  (c2-C;i)T)Vt(C  -  (c2-c1)T) 

=  H  (T)V  (£) 

1 

so 

V(£>  =  Hc  (T)Vt(C  -  (c2-Cl)T) 

this  implies  v(£)  e  w  (U  u  V  ) ,  but  by  the  same  argument  as  above 

C2  °1  C1 

uj  (u  u  V  )  =  u  (U  u  V  )  =  w  (U  u  V)  .  So  v(£)  e  D  implies 
C2  C1  C1  °2  °2  C2  °2  C1 

v(£)  e  D  ifc^ci=  1,2,  call  this  common  set  D.  We  must 
C2 

show  that  D-  =  D .  Let  u(C,t)  be  an  orbit  in  D_,  then  either 
c  c 

u(£,t  )  e  U_  n  V_  for  a  sequence  t  u(£,t)  e  V  for  all  large 

n  c  c  n 

negative  t  or  u(£,t)  e  V  for  all  large  negative  t.  In  the  former 

case,  it  is  a  simple  consequence  of  Theorem  5.1  that  u(£,t)  would  be 

the  travelling  wave  and  so  in  D.  In  either  of  the  latter  cases 

H  (t)u(£,0)  would  lie  in  CJ  or  V  respectively  for  some  c  /  c  and 
c 

large  negative  t.  This  is  either  impossible  or  it  puts  u(C»t)  in  D. 

The  proof  will  be  complete  if  we  can  show  that  D  =  W.  We  con¬ 
sider  D  in  the  c-semiflow  (for  the  sake  of  notation,  we  shall  drop 
the  c  in 

c 


-  of*.-- 


Define  the  set  s(U),  for  a  given  set  U,  to  be  the  set  of  u('l 

where  u(£+s)  e  u.  We  firstly  show  that  if  u  e  D  either  (1) 

u  c  s(Uc  n  V  )  ,  for  some  s,  or  (2)  u(t,)  <_  e  for  all  "  or  (3) 

u(£)  >_  6  for  all  £ . 

Suppose  neither  (2)  nor  (3)  are  true  then  there  is  a  '  so  that 
e  <  u(C)  <  <5.  It  must  be  true  that  lim  u(t)  <  t  and  lim  u(',)  •, 

otherwise  some  translate  of  u(£)  would  not  lie  in  U  i  V  and  this 

c  c 

is  impossible  as  D  is  translation  invariant  (argued  above  for  , 

c  ^  c)  and  D  c  U-  l  V-  . 

c  c 

It  follows  that  the  set 

G  =  {p  •.  33  |  u(£)  c  s  (V-)  } 

is  non-empty.  Let  s  =  inf  G,  then  s  •  G  as  V  is  closed.  We 

c 

claim  that  u(£)  c  s(U-),  if  this  were  not  true  there  would  exist  a 

c 

p  so  that  both  u(£)  j  p(V-)  and  u(C)  p(U-)  by  minimality  of  s 

c  c 

and  the  fact  that  U-  is  closed.  But  this  is  a  contradiction  as  we 

c 

would  then  have  u(£+p)  i  U-  u  V- ,  and  u(£+p)  £  D  c  U-  u  V- . 

c  c  c  c 

So  one  of  the  alternatives  (1),  (2)  or  (3)  above  is  true  of 

u(£)  e  D.  If  either  (2)  or  (3)  were  true,  since  D  is  invariant,  we 

would  have  u(£)  e  {0,1}.  So  we  have  shown  that 

(5.6)  D  c  {0,1}  u  (  U  s(U-  n  V-)) 

_  c  c 
se  33 

But  it  follows  easily  from  theorem  5.1  that  w-(  U  s(U-  n  V-))  -•  W 

c  _  c  c 
se  33 

and  so  w-(D)  =  W.  Since  w-(D)  =  D,  D  =  W. 
c  c 


0S..-W  • 


★ 

From  above  D  =  ^  (U  i  V)  for  anv  c  -  i  3 ,  c 
c 

the  proof  of  the  theorem. 


Remark:  (1)  In  lemmas  5.2  and  5.3,  0  and  1  are  r.c t  n  :  •  .  ■ 

for  the  case  c  =  c,  this  is  because  they  obviously  ar<  -  :  •  .  • 

in  that  semiflow,  the  travelling  wave  is  a  string  of  _-r  i"  i  -a:  ■  •  • 

that  approaches  both  of  them.  (2)  From  Lemmas  5.1,  3.2  and  T:.\  :■  :• 
we  have  all  we  need  for  the  theorems  of  chapter  4,  section  III 
specialising  to  the  case  c  =  0.  Mote  however  that  we  did  nor  — 
this  case  in  isolation,  the  above  proof  depended  crucially  on  us.:-:  •: 
spectrum  of  semiflows. 


III.  APPLICATION  OF  EXTENDED  CONTINUATION 


Again  we  consider  the  bistable  equation  (1.3)  where  f  satis  :i-. 
(HI)  ( (H2)  is  not  needed  in  this  section  at  all). 

We  will  apply  the  semiflow  of  section  I  and  the  results  of  secti  r. 
II  to  the  question  of  how  solutions  behave  when  followed  out  in  a 
radial  direction  with  speed  c.  If  u(r,t)  is  a  solution,  for  a  <;i ••••■:. 
c,  we  want  to  study  its  behaviour  along  lines  r  -  ct  =  constant, 
for  any  p  e  B  set  r  =  p+ct  and  we  get  a  function  v(p,t)  = 
u(p+ct,t)  whose  domain  is  time  dependent  i.e.  p  >_  -ct.  However,  for 
each  p  e  B  there  is  a  T  so  that  if  t  >_  T  v(p,t)  is  defined  and; 
so  it  makes  sense  to  try  and  determine 


i- 


* 


i U ( a , e )  l  V(b,e) ]  *  [0,c  ] 

where  a,  c,  b,  6  satisfy  all  the  conditions  given  in  section  II. 

We  will  compute  (5.7)  for  initial  data  u(r)  •  A^  that  satisfy 

(5.8)  lim  u(r)  <  a 
r^+co 

(5.9)  u(r,t)  ->-1  as  t  -»  °°,  uniformly  on 
compact  sets. 

Theorem  5.3.  Suppose  u(r)  e  satisfies  (5.8)  and  (5.9)  then  there 

is  a  function  c(t)  =  o(t)  so  that 

(5.10)  lim  ||u(p+ct+  ~(t)  ,t)  -  u  (p)  ;|  ^  =  0 

where  u^  is  the  one-dimensional  travelling  wave  and  its  speed  is  c. 

Remarks:  (1)  The  theorem  says  that  any  spherically  symmetric 

solution  (with  the  restriction  of  being  in  AQ)  which  propagates 
(condition  (5.6))  takes  the  shape  of  the  one-dimer.sional  travelling 
wave  as  t  -*■  +°°. 

(2)  The  function  u (p+ct+ v( t) , t)  is  actually  defined  on 
a  time-dependent  domain  of  p's.  For  the  norm  to  make  strict  sense  we 
can  extend  the  function  to  be  constant  for  p  <_  -ct-s:(t)  . 

(3)  An  alternative  way  to  stating  (5.10)  would  be 

(5.11)  lim  ||  u  (r ,  t)  -  u  (r-ct- v?(t)  )  ||  =  0  , 

■I  I  1  oo 

t-*» 

where  the  supremum  is  now  taken  over  r  e  3R+  .  This  is  how  the  state¬ 
ment  of  theorem  5.1  was  made.  It  is  clear  that  (5.10)  and  (5.11)  are 


equivalent . 


1 
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The  main  idea  in  the  proof  of  theorem  5.3  is  a  shooting  argument 

for  the  semiflow  H(t).  From  knowing  the  behaviour  of  lim  H  (t)u 

t-H-00 

for  c  <  c  and  c  >  c  we  shall  make  a  conclusion  about  lim  H-(t)u. 

c 

t-*-+°° 

(Notation:  we  shall  denote  an  element  of  A  by  the  component  that  is 

a  function  in  M,  the  X -co-ordinate  is  implicit  in  the  symmetry  of 
this  function.) 

We  will  need  a  perturbation  statement  similar  to  theorem  4.1  but 
for  this  extended  semiflow  on  A  x  [0,“).  A  neighbourhood  of  a  function 
v(x)  €  A^,  depending  on  a  compact  set  K,  contained  in  the  x^-axis, 
and  an  £  >  0,  is  of  the  form 

(5.12)  {u(x)  |  |u(x)-v(x) |  £  e  for  x  e  K} 

Let  U  be  such  a  set,  or  a  union  of  such  sets.  Let  U  c  M  be  defined 
by  the  same  inequalities.  Set 

(5.13)  U  =  (U  x  [0,1])  n  A 

(5.14)  UX  =  U  n  (Mx  [ X , 1] )  . 

(It  is  to  be  understood  that  if  a  Roman  letter  refers  to  a  set,  the 
corresponding  script  letter  refers  to  the  associated  set  given  by  (5.13). 
Similarly,  if  a  superscript  X  appears.) 

Let  I  c  [0,<”)  be  compact  and  N  =  Uxicax  [0,°°)  .  N  and 
have  their  obvious  meanings.  In  the  following  refers  to  the 


action  of 

H  ( t) 

on  A  x  I . 

Lemma 

5.4. 

If,  for  some 

x  >  0 

(5.15) 

cl  (N  • 

[t,°°))  c  Int(N) 

then  there  is  a  X  <  1  so  that 


(5.16)  ci.(NX  •  [t,*)1  c  Int  (N ' )  . 

Proof  of  Lemma  5.4.  It  suffices  to  show  that  for  each  , 

there  is  a  X  <  1  so  that 

(5.17)  cK(NX  •  [t  ,T  ])  ■=  Int  (!)')  . 

Setting  =  t,  =  2t  in  (5.17)  and  iterating  yields  (5.16)  . 

We  can  pick  an  open  set  V  and  a  closed  set  C,  in  A^,  so  that 
V  c  c  c  int(N)  and 

(5.18)  cl (N  •  [t  ,T  ])  c  v  . 

Now  suppose  there  is  no  X  for  which  (5.17)  is  satisfied.  Then 

X 

there  is  a  sequence  {u  }  with  u  e  U  ,  X  -*■  1,  sequences  {c  - 

n  n  n  n 

and  {t  }  so  that  (u  ,c  )  •  t  i  K. 
n  n  n  n  ^ 

H(t)  is  a  compact  semiflow  as  compactness  was  proved  by  a  global 

derivative  estimate  which  is  obviously  preserved  under  translation.  It 

follows  that  { (u  ,c  )  •  t  }  has  a  limit  point.  We  can  assume  that 
n  n  n 

(u  ,c  )  •  t  ■*  u,  c  -*•  c*  and  t  •+  t'.  Using  the  fact  that 
n  n  n  n  n 

{  (un'cn)  *  is  precompact,  it  is  a  standard  argument  to  see  that 

(u  ,c  )  •  t  and  (u  ,c‘)  •  t’  have  the  same  limit, 
n  n  n  n 

Let  {v  }  be  a  sequence  in  U  so  that  u  -v  ->-0,  for  instance 

n  n  n 

let  v  agree  with  u  on  the  x-axis  above  -X  /1-X  .  In  theorem 

n  n  1  n  n 

3.2  we  actually  proved  that  the  semiflow  is  uniformly  continuous  with 
respect  to  initial  data,  this  is  a  consequence  of  estimate  (3.16). 


It  follows  that  this  is  also  true  for  H(t)  on  A. 


It  is  then 


check  that  (v  ,c')  •  t*  -  u,  but  u  i  V  and  so  this  contradicts 

n 

(5.18)  . 

Proof  of  Theorem  5.3.  Let  u(r)  e  satisfy  (5.8)  and  (5.9).  Kecai;- 

0 

*  — 

ing  c  from  section  II,  if  is  between  0  and  c,  set 

* 

I  =  [c^c  ]  and  E  =  u  *  I. 

There  is  a  set  N  =  (U  t  V)  *  I  which  is  an  attracting  neighbour¬ 
hood  for  the  two  dimensional  invariant  manifold  of  Figure  5.4.  By 

Lemma  5.4,  there  is  a  X  <  1  and  x  >  0  so  that 

X  x  *  X 

c£(N  •  [x,”))  c  Int(N).  Since  c  >0,  u(N  )  c  A^  and  so  N  is  an 

attracting  neighbourhood. 

If  b  is  large  enough  T (b) E  c  by  condition  (5.8)  .  Conse- 

* 

quently  w(T(b)E)  c  w  x  [ct,c  ]. 

Let  d  be  a  metric  on  M.  Pick  any  6  so  that  0  <  6  <  d(0,l) . 

The  sphere  S  =  {u  |  d(u,l)  =  6}  separates  M  into  two  open  sets,  its 

* 

interior  and  exterior  D  =  (S  x  [0,1])  x  [c*,c  ]  also  separates 

* 

A  x  [c^,c  ]  into  two  disjoint  open  sets,  with  0  in  one  of  them  and 
1  in  the  other.  If  we  can  show  that  for  large  t  H(t)(T(b)E)  inter¬ 
sects  both  these  sets,  then  by  connectedness  it  must  intersect  D. 

We  will  show  that  for  any  fixed  b  0 


(5.19) 


fl 

lim  H  (t)  (T(b)  u)  -l 
t-H-<»  C  I 

lo 


if  c  <  c 


if  c  >  c 


If  (5.19)  is  true  for  one  particular  b  it  is  true  for  every  b,  by 
translation  invariance.  To  prove  the  second  part,  with  the  c-semiflow 

l 

I 


on  A,  0  is  an  attractor  by  lemma  5.1  with  U  as  an  attracting 
neighbourhood.  By  lemma  5.4  there  is  a  U  which  is  an  attracting 
neighbourhood  for  0.  If  b  is  large  enouah  T(b)u  r  U  ana  so  the 
limiting  behaviour  follows. 

For  the  first  part,  consider  the  set  G  =  (u(r,t)  ;  u(r,0)  =  u(r) 

By  lemmas  5.2  and  5.4  there  is  an  attracting  neighbourhood  of  1,  V', 
for  the  c-semiflow  if  0  <  c  <  c.  If  b  is  large  enough 
T(b)G  n  V  /  i ,  since  1  e  c£(T(b)G)  by  (5.9)  and  translation  invari¬ 
ance.  But  then  lim  H  (t) (T(b)g)  =  1  for  some  g  €  G.  Since 

c 

t-H-00 

g  =  u(r,t),  for  some  T,  the  first  part  of  (5.19)  must  hold. 

We  now  know  that  for  each  large  •«b,  *f("t)' (Ttbf£)’  '  intersects  D. 

Then  there  is  a  function  c(t)  such  that  H  ,^(T(b)u)  e  D.  As 

c  (t) 

★ 

t-M-oi,  h  (T(b)u)  -*■  W  x  (c  ,c  ]  ,  but  D  is  closed  and  there  is  only 
c  ( t )  * 

one  point  in  W  that  is  a  distance  6  from  1,  so  H  ,  , (T(b)u)  con- 

c  (t) 

verges  to  this  point,  which  is  not  0  or  1  and  therefore  is  a 
travelling  wave,  call  it  u  . 

Rewriting  this  quantity  by  restricting  to  the  x^-axis  and  letting 
y (t)  =  c (t) t,  we  get 

(5.20)  lim  u(x1  +  y(t),t)  =  u^x^) 

t-H-® 

(5.20)  is  actually  in  the  compact-open  topology  but  since  everything  is 
nonincreasing  in  x^  (if  it  is  >_  -y(t))  then  (5.20)  also  holds  in 
the  sup-norm.  Letting  y(t)  =  ct  +  *fi( t)  and  p  =  x^  we  have  (5.10) . 
That  <?(t)  =  o(t)  is  an  easy  consequence  of  (5.19). 


Remarks:  (1)  What  is  most  interesting  in  Theorem  5.3  is  that  or.’, 

knowledge  of  one-dimensional  behaviour  is  used.  We  prove  nothir.cs 
directly  about  the  spherically  symmetric  solutions  except  by  perturba¬ 
tion  from  their  one-dimensional  limiting  behaviour. 

(2)  It  is  not  hard  to  see  that  c(t)  -*  -°°  as  t— *-»,  so  the 
spherical  wave  may  lag  behind  the  one-dimensional  wave,  but  does  not 
get  ahead.  This  follows  from  the  fact  that  any  solution  to  the  one¬ 
dimensional  equation  is  a  supersolution  for  the  spherically  symmetric- 
equation  and  such  a  solution  can  easily  be  found  which  majorises  u(r) 
and  tends  to  the  travelling  wave  (by  Theorem  5.1) . 


Chapter  6 


Conclusion 

I.  INTERPRETATION  OF  RESULTS 

In  the  first  two  sections  of  this  chapter,  we  shall  discus- 
sively  the  bistable  equation  (1.3)  with  the  nonlinearity  satisfy::. 
(HI)  and  (H2 ) . 

At  the  end  of  chapter  4  we  found  a  set,  called  VJ  ,  that  is  c. 
attractor  in  the  semiflow  on  A^  (the  world  of  spherically  symmetric 
data)  .  It  is  perhaps  not  clear  what  the  fact  that  this  set  is  ar. 
attractor  means  for  the  behaviour  of  the  solutions  of  the  equation, 
this  section  we  shall  give  an  interpretation  of  this  statement  in  r.or 
traditional  terms. 

Recall  that  AQ  consists  of  functions  in  B  that  are  spherical 
symmetric,  between  0  and  1,  nonincreasing  in  r.  Define  the  set 

U  =  (u(r)  e  Aq  |  u{+®)  <  a} 

being  an  attractor  translates  into  the  following  theorem. 

Theorem  6.1.  U=U_uU-uU_  where 
-  0  u  1 


=  tu(r) 

|  u(r,t) 

0 

as 

t  00 

uniformly) 

=  (u(r) 

|  u(r,t) 

-*■  u 

as 

t  00 

uniformly } 

=  {u (r) 

|  u(r,t) 

-*•  1 

as 

t  -►  » 

uniformly 

on  compact  sets} 


Further  Uq  and  are  both  open  and  connected. 
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Proof.  Given  any  u(r)  U,  there  is  an  a  and  an  z  so  that 
u  c  U  (a,e),  pick  a  corresponding  V  (b,5)  then  U  i_  V  is  an 

5  S  S  s 

attracting  neighbourhood  for  W 


and  so  everything  in  U  'J  tends  to  one  of  0,  u  and  1/  there¬ 

fore  so  does  u(r) ,  it  is  obvious  that  in  the  first  two  cases  the  con¬ 
vergence  is  uniform. 

Since  0  and  1  are  attractors  in  their  own  right,  rJQ  and  U 
must  both  be  open  and  since  neighbourhoods  of  0  and  1  are  connected 
they  must  also  be  connected.  This  completes  the  proof. 

Since  U  is  obviously  connected  the  t'neoren  says  that  the  removal 

of  U-  disconnects  the  set  into  these  two  components  that  give  the 
u 

propagation  and  decay  behaviour  respectively.  So  separates  the 

two  regimes  and,  in  some  sense,  gives  the  optimal  set  of  comparison 
functions  for  both  propagation  and  decay.  This  is,  practically  speak¬ 
ing,  not  particularly  useful  as  these  are  hard  to  locate  except  near  u 
but  they  do  give  an  exact  picture  of  the  situation. 


Another  way  to  express  these  facts  is  that  U-  is  a  codinension 
one  set  and  so  it  is  "small".  The  two  pieces  it  separates  are  the 
regimes  of  propagation  and  decay  respectively. 

II.  ASYMPTOTIC  STATES  AND  PERMANENT  SOLUTIONS 

We  have  not  quite  shown  that  the  propagation  regime  U^  is  an  SAS 
with  respect  to  A^  (see  definition  1.4).  It  is  open  in  the  compact- 
open  topology  and  so  also  in  the  sup-norm  but  it  would  require  an  esti¬ 
mate  on  how  ^(t)  varied  across  initial  data  to  show  that  elements  of 
had  the  same  sup-norm  asymptotic  behaviour,  (5.10)  is  not  quite 
enough.  If  a  definition  were  formulated  which  bore  the  relation  to 
definition  1.2  that  orbital  stability  does  to  stability  in  qualitative 
ordinary  differential  equations,  would  be  an  SAS  in  this  other  sense. 

Aq  has  three,  albeit  invariant,  conditions  on  it.  As  far  as  U^ 
being  an  asymptotic  state  is  concerned,  it  is  not  hard  to  drop  the  con¬ 
dition  0  _<  u(r)  <_  1,  and  it  may  be  possible  to  drop  the  monotonicity 
condition.  But  the  major  open  question  is  whether  the  condition  of 
spherical  symmetry  is  necessary.  For  instance  it  is  not  known  whether 
there  could  be  a  wave  that  propagates  in  IR2  with  an  elliptical 
wave-front. 

For  the  above,  the  only  feature  of  used  is  that  it  contains 

three  equilibrium  points  and  orbits  joining  them.  We  did  not  use  the 
fact  that  there  are  only  two  connecting  orbits  in  it.  This  was  proved 
in  chapter  4.  In  terms  of  theorem  6.1,  there  is  one  in  UQ  and  one  in 


U^.  Each  of  these  connecting  orbits  is  a  permanent  solution  (see  defi 
nition  1.3)  .  Consequently  there  is  a  unique  permanent  solution  in  oac 
SAS  UQ  and  . 

These  permanent  solutions  are  "special"  solutions  of  the  equation 
as  we  would  generally  only  expect  forward  existence  in  parabolic  equa¬ 
tions.  Nevertheless,  it  is  not  clear  what  the  real  significance  of 
their  being  unique  is.  Such  permanent  solutions  as  equilibria  and 
travelling  waves  supply  the  asymptotic  information  for  their  SAS.  If 
we  consider  the  SAS  this  information  is  given  by  the  one-dimen¬ 

sional  travelling  waves,  not  this  permanent  solution  in  U  .  This 
solution  may  however  be  relevant  to  the  asymptotic  behaviour  of  non- 
spherically  symmetric  data  since  it  contains  spherical  information 
and  the  travelling  waves  do  not. 

III.  COMMENTS  ON  THE  TECHNIQUE 

The  main  technique  in  this  work  is  the  construction  of  the  space 
A  which  relates  the  one-dimensional  and  spherically  symmetric 
behaviour  of  a  reaction  -  diffusion  equation.  The  space  realises  the 
intuition  of  the  one-dimensional  world  living  out  "at  infinity"  in  the 
spherical  world.  By  constructing  a  concrete  mathematical  object,  such 
as  A,  we  can  see  exactly  what  is  forced  on  the  spherical ly-symmetric 
solutions  by  this  relationship. 

The  presentation  of  this  construction  here  treats  this  as  a  homo- 
topy-continuation  method.  But,  in  many  respects,  it  is  more  closely 
related  to  the  technique  in  'Ordinary  Differential  Equations'  of 


pasting  a  manifold  onto  phase  space  which  represents  son-  I:r. :  *;• 
behaviour.  From  studying  the  induced  flow  on  this  rnanifci  i,  : 

sions  can  be  made,  by  continuity,  about  the  flow  on  thv  or :  a •  i 
space.  A  recent  example  of  the  power  of  this  technique  is  1:. 

[1],  and  this  work  is  partially  responsible  for  ir.st-ir mcs  t 
this  idea  here. 

It  is  a  drawback  in  the  technique  that  the  underlying  tip  1  : 

A  is  the  compact-open  topology  because  stability  statements  are  ~ 
easily  made  in  the  sup-norm  topology.  This  topology  is  used  in  rv.v 
essential  ways,  for  compactness  of  the  semiflow  and  to  make  the  spare 
tie  together  so  that  perturbation  is  possible  from  A^  to  A. . 

The  way  to  circumvent  the  problem  is  to  restrict  the  semiflow  to 
an  invariant  subspace.  This  subspace  should  have  the  property  that  an 
open  neighbourhood  of  the  solution  whose  stability  is  in  question  with 
the  inherited  topology  should  be  effectively  the  same  as  that  set  with 
the  sup-norm  topology.  For  example,  if  the  stable  one-dimensional 
behaviour  expected  is  given  by  a  travelling  pulse  or  front  then  the 
invariant  subspace  should  consist  of  functions  that  look  like  the  wave 
as  x  -+  ±°°.  So  the  invariant  subspace  makes  the  data  be  close  to  the 
wave  at  ±°°  and  consequently  being  close  in  the  compact-open  topoloqy 
implies  the  same  in  the  sup-norm  topology.  This  also  shows  the 
necessity  of  considering  an  interval  of  different  speeds.  The  wave  car* 


move,  metaphorically,  "out  of  sight" 


in  the  compact-open  topoloqy  units. 


we  focus  on  it  at  the  correct  speed. 


i 


So  the  key  to  a; -plying  the  technique  lies  ir. 
right  places  where  the  space  should  be  tightened 
logy  or  loosened  to  the  c omp a c t - o ;>e n  topology. 
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